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I.  INTRODUCTION 

A.  PURPOSE  AND  OVERVIEW 

The  main  purpose  of  this  thesis  is  to  conduct  an  experiment,  using  two  different 
types  of  combat  attrition  models,  to  determine  whether  their  results  are  similar.  There 
is  a  belief  that  deterministic  Lanchester  Equation  type  models  closely  parallel  the  mean 
results  of  stochastic  computer  simulations.  This  research  will  attempt  to  determine 
whether  the  results  of  a  deterministic  "Generalized  Lanchester  Equation"  type  model 
and  a  stochastic  computer  simulation  model  are  similar.  This  work  is  the  first  such 
research  of  this  type  published  to  the  author's  knowledge. 

Chapter  One  provides  the  reader  with  an  introduction  to  modeling  and  an 
overview  of  the  history  of  Lanchester  Equation  type  models.  The  purpose  of  Chapter 
Two  is  to  introduce  the  two  models  which  are  utilized  in  this  experiment  and  Chapter 
Three  is  a  presentation  of  the  experimental  design.  The  results  of  the  experiment  are 
provided  in  Chapter  Four  and  a  discussion  of  the  conclusions  and  recommended 
follow-up  research  are  given  in  Chapter  Five. 

B.  MODELING 

Throughout  history,  military  analysts  have  been  developing  and  utilizing  various 
types  of  models  for  detailed  analysis  to  assist  in  decision  making.  A  model  is  defined 
as  a  "simplified  representation  of  the  entity  it  imitates  or  simulates"  [Ref.  1:  p.  I].  A 
second  definition  of  a  model  used  by  the  U.S.  Army  Models  Review  Committee,  "an 
abstract  representation  of  reality  which  is  used  for  the  purpose  of  prediction  and  to 
develop  understanding  about  the  real-world  processes"  [Ref.  2:  p.  5],  implies  that 
models  are  designed  to  be  close  representations  of  a  real-world  entity  or  process.  A 
final  definition  for  a  military  model  is,  "an  abstraction  of  reality,  the  elements  of  which 
are  chosen  for  (a)  an  investigative  purpose  or  (b)  a  resource  management  purpose;  in 
other  words,  an  abstraction  to  assist  in  making  decisions."  [Ref.  1:  p.  3].  From  these 
definitions,  it  can  be  seen  how  important  the  development  and  use  of  models  is  to  the 
military  analyst. 

In  general,  there  are  three  types  of  models: 

1)  Iconic  model 

2)  Analogue  model 

3)  Symbolic  model 

12 


An  iconic  model  is  a  miniature  version  of  the  entity,  such  as  an  airplane  or  a  tank 
model.  An  analogue  model  is  an  artificial  representation  of  reality.  An  example  of  an 
analogue  model  is  a  map  which  represents  the  three  dimensional  real-world  on  a  two 
dimensional,  small  scale  sheet  of  paper.  A  symbolic  model  is  one  in  which  words  or 
numerical  descriptions  are  used  to  represent  an  entity  or  process.  An  example  of  a 
symbolic  model  is  a  mathematical  equation  which  represents  a  process  such  as 
attrition. 

The  types  of  models  used  by  military  analysts  cover  the  entire  spectrum  from 
military  field  exercises  (which  can  be  thought  of  as  iconic  models)  to  analytical  models. 
Figure  1.1  is  a  combination  of  information  found  in  Taylor's  book  Force-on-Force 
Attrition  Modelling  [Ref.  2:  p.  7],  and  Hughes'  book  Military  Modeling  [Ref.  1:  p.  10], 
which  shows  the  various  types  of  models  in  use  and  their  characteristics  of  operational 
realism,  degree  of  abstraction,  convenience  and  accessibility. 


MILITARY  MILITARY  MAP  WAR  COMPUTER    ANALYTICAL 

FIELD  FIELD  EXERCISES        GAMES      SIMULATIONS      MODELS 

EXERCISES  EXPERIMENTS 


HUMAN    DECISION    IMPACT   AND    OPERATIONAL    REALISM   \ 


DEGREE    OF    ABSTRACTION 


OUTCOME     REPRODUCIBILITY 


CONVENIENCE    AND    FLEXIBILITY 


Figure  1.1     Current  Model  Types. 
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The  three  most  convenient  and  accessible  types  of  models  shown  in  Figure  1.1 
are  war  games,  computer  simulations  and  analytical  models.  These  three  types  of 
models  are  also  the  most  commonly  used  to  model  combat  attrition.  A  brief 
description  of  each  of  them  will  be  given  to  better  understand  and  help  distinguish 
between  them. 

1.  War  Games 

War  games  are  conducted  using  real-world  combat  scenarios  and  people 
playing  various  positions  on  a  headquarters  battle  staff.  A  war  game  is  conducted  to 
allow  a  commander  and  his  staff  to  exercise  the  staff  decision  making  processes  in 
'realistic'  combat  scenarios.  This  allows  the  commander  to  exercise  and  develop  his 
staff  without  incurring  the  true  cost  of  combat,  the  loss  of  human  life.  War  games  can 
be  run  utilizing  the  sand  table  method  or  on  computers.  The  key  to  war  gaming  is  that 
there  is  always  a  person  making  the  decisions  that  are  required  in  the  game. 

2.  Computer  Simulations 

Computer  simulations  are  also  used  to  model  the  combat  attrition  process. 
The  primary  difference  between  computer  simulations  and  war  games  is  the  method 
used  to  represent  the  decision  process.  War  games  use  people  playing  staff  roles  to 
make  the  decisions  where  simulations  use  algorithms  to  represent  the  decision  process. 
A  combat  simulation  not  only  represents  the  combat  process,  it  acts  it  out  from  start 
to  finish.  A  computer  simulation  begins  with  a  set  of  input  parameters  and  runs 
continuously  until  completion  of  the  battle  and  then  provides  the  results  of  the  battle 
as  output.  Computer  simulations  which  utilize  pseudo-random  number  generators  to 
determine  the  results  of  random  events,  such  as  the  outcome  of  one  soldier  firing  at 
another,  are  called  Monte  Carlo  Simulations  and  are  stochastic  models. 

3.  Analytical  Models 

Analytical  models  are  the  third  type  of  model  used  in  modeling  combat 
attrition.  Analytical  models  are  symbolic  models  which  use  mathematical  symbols  and 
equations  to  represent  the  combat  attrition  process.  Analytical  models  can  be 
developed  as  stochastic  models  or  deterministic  models.  Stochastic  models,  as 
introduced  in  the  previous  section,  utilize  probablity  distributions  to  determine  certain 
variables.  Therefore  the  output  from  two  consecutive  runs  of  a  stochastic  model  will 
more  than  likely  be  different.  The  deterministic  model  however,  will  produce  the  same 
set  of  output  values  for  a  given  set  of  input  parameters. 
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My  research  is  focusing  on  the  comparison  of  a  stochastic  Monte  Carlo 
Simulation  and  a  deterministic  analytical  model.  The  deterministic  model  I  am  using 
for  comparison  is  a  Generalized  Lanchester  Equation  model  which  was  developed  from 
the  original  Lanchester  Equation  Theory.  To  provide  a  basis  for  the  readers 
understanding  of  this  experiment,  I  will  give  a  brief  introduction  to  the  history  of  the 
Lanchester  Equation  type  models. 

C.   LANCHESTER  EQUATION  THEORY 

1.  History  Of  Lanchester  Equation  Type  Modeling 

F.  W.  Lanchester  was  a  British  aeronautical  engineer,  who  in  1914  developed 
two  sets  of  simple  differential  equations  to  model  the  combat  attrition  process  between 
two  opposing  homogeneous  forces  (i.e.,  fighter  aircraft  vs..  fighter  aircraft).  His 
intentions  were  to  provide  insight  into  the  dynamics  of  combat  under  'modern 
conditions'  of  warfare  and  to  justify  the  principle  of  concentration  of  forces.  The  term 
modern  conditions  of  warfare  will  be  addressed  again  later. 

Lanchester's  original  work  was  designed  to  model  a  force-on-force  attrition 
process  involving  two  homogeneous  forces.  He  developed  two  combat  attrition 
models,  each  having  it's  own  unique  assumptions,  for  this  purpose.  These  original 
models  were  designed  to  model  an  'aimed  fire'  combat  scenario  and  an  'area  fire' 
combat  scenario.  These  two  models  are  often  referred  to  as  the  classical  Lanchester 
Equation  Theory. 

2.  Modeling  Aimed  Fire 

The  phrase  aimed  fire,  as  Lanchester  used  it  in  his  original  work,  refers  to 
combat  between  two  forces,  X  and  Y,  where  each  combatant  from  X  force  acquires  a 
Y  force  target  (i.e.,  locates  and  takes  aim)  and  fires.  An  example  of  this  type  combat  is 
an  infantry  battle  for  control  of  an  area  such  as  hill  224.  The  following  assumptions 
apply  for  this  model: 

1)  The  combat  being  modeled  involves  two  homogeneous  forces  (i.e.,  Infantry  vs. 
Infantry). 

2)  The  entire  X  force  and  Y  force  are  within  weapons  range  of  one  another. 

3)  The  effects  of  weapons  rounds  are  independent. 

4)  Each  of  the  forces  is  well  enough  aware  of  the  location  and  condition  of  all 
enemy  forces  so  that  they  will  engage  only  live  enemy  units.  Also  the  rate  at 
which  they  kill  enemy  targets  is  constant. 

5)  Fire  is  uniformly  distributed  over  surviving  enemy  targets. 
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Lanchester's  original  model  for  aimed  fire  is  given  by  Equation  Set  1.1: 

dx/dt  =  -ay(t)  withx(0)=  xQ, 

(eqn  1.1) 
dy/dt  =  -bx(t)  withy(0)=  yQ, 

where  x(t)  is  defined  as  the  number  of  combatants  of  X  force  at  time  t,  y(t)  is  defined 
as  the  number  of  combatants  of  Y  force  at  time  t,  and  a  and  b  are  constants  which  are 
called  attrition-rate  coefficients.  The  attrition-rate  coefficient  'a'  represents  the 
effectiveness  of  a  Y  combatant  killing  a  X  combatant  per  unit  of  time.  Similarly  the 
attrition-rate  coefficient  'b'  represents  the  effectiveness  of  a  X  combatant  killing  a  Y 
combatant  per  unit  of  time.  This  set  of  equations  assumes  that  the  rate  at  which  a 
force  X  can  cause  casualties  to  a  Y  force  is  proportional  to  the  number  of  combatants 
in  the  X  force.  Similarly,  the  rate  at  which  a  combatant  from  Y  force  can  cause 
casualties  to  the  X  force  is  assumed  to  be  proportional  to  the  number  of  combatants  in 
the  Y  force.  When  Equation  Set  1.1  is  integrated,  the  resulting  state  equation  is 
obtained  and  has  been  labeled  the  Lanchester  Square  Law: 

b(x02-x(t)2)  =  a(y02-y(t)2).  (eqn  1.2) 

3.  Modeling  Area  Fire 

The  phrase  area  fire  as  Lanchester  used  it  in  his  original  work,  refers  to  a 
combat  scenario  involving  two  homogeneous  forces,  say  X  and  Y,  which  are  uncertain 
of  exact  enemy  locations  and  engage  one  another  by  firing  in  the  general  area  where 
the  enemy  force  is  located.   This  scenario  has  the  following  assumption  set: 

1)  The  combat  being  modeled  involves  two  homogeneous  forces  (i.e.,  Infantry  vs. 
Infantry) 

2)  The  entire  X  force  and  Y  force  are  within  weapons  range  of  one  another. 

3)  The  effects  of  weapons  rounds  are  independent. 

4)  Each  force  is  aware  only  of  the  general  location  of  the  enemy  force  and 
therefore  engages  the  enemy  by  firing  into  that  general  'area'  without  the 
benefit  of  knowing  their  effectiveness. 

5)  Fire  from  all  surviving  combatants  is  uniformly  distributed  over  the  area 
which  the  enemy  occupies. 

6)  Each  force  has  the  same  vulnerable  area  to  enemy  fire. 
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Lanchesters  original  model  for  area  fire  is  defined  by  the  Equation  Set  1.3: 

dx/dt  =  -axy, 

(eqn  1.3) 
dy/dt  =  -byx. 

This  model  assumes  that  the  rate  of  attrition  is  proportional  to  both  force  levels  x(t) 
and  y(t).  This  means  that  in  the  area  fire  model,  the  rate  at  which  a  X  force  can  cause 
casualties  to  a  Y  force  is  not  only  dependant  upon  the  number  of  combatants  in  its 
own  force,  but  also  on  the  number  of  combatants  alive  in  the  Y  force.  The  rate  at 
which  the  Y  force  can  cause  casualties  on  the  X  force  is  similarly  dependent  on  both 
their  own  force  level  and  the  force  level  of  the  opposing  force  X. 

The  well  known  'Linear  Law"  is  obtained  by  integration  of  Equation  Set  1.3 
and  is  shown  here  in  Equation  1.4  : 

b(x0-x(t))=a(y0-y(t)),  (eqn  1.4) 

where  xQ  and  y0  are  the  initial  force  levels  at  time  t  =  0,  x(t)  and  y(t)  are  the  force 
levels  at  time  t,  and  a  and  b  are  the  attrition-rate  coefficients  for  each  force. 
4.  Attrition  Rate  Coefficients 

The  Lanchester  attrition-rate  coefficients  were  introduced  in  the  models  above 
and  were  defined  as  constants.  In  order  to  use  the  Lanchester  models  there  must  exist 
a  set  of  acceptable  attrition-rate  coefficients  for  each  possible  type  of  engagement. 
Lanchester's  original  model  utilizing  homogeneous  forces,  only  requires  the  availability 
of  a  limited  number  of  coefficients.  However,  in  Lanchester  Equation  type  models 
which  allow  for  the  combat  between  heterogeneous  forces,  the  number  of  attrition-rate 
coefficients  grows  very  rapidly  when  one  considers  all  the  possible  combinations  of 
combatant/weapons  vs.  combatant/weapons  that  can  occur.  As  Lanchester  postulated 
his  original  work,  he  used  constants  for  attrition-rate  coefficients,  which  implies  that 
the  kill  capability  of  a  force  does  not  change  over  time.  With  constant  attrition-rate 
coefficients  these  models  are  easy  to  solve.  However,  there  is  no  reason  why  we  should 
believe  that  such  factors  as  range,  weather,  visibility,  and  training  do  not  affect  the  kill 
capability  of  a  combatant.  Taylor  provides  two  methods  for  determining  numerical 
values  for  Lanchester  attrition-rate  coefficients  which  are  used  in  the  United  States. 
[Ref.  2:  p.  45] 
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1)  A  statistical  estimate  based  on  'combat'  data  generated  by  a  detailed  Monte 
Carlo  combat  simulation. 

2)  An  analytical  submodel  of  the  attrition  process  for  the  particular  combination 
of  firer  and  targets. 

In  the  first  method,  the  simulation  output  is  used  to  fit  one  or  more  free 
parameters  in  the  analytical  model  in  an  attempt  to  have  it  to  provide  similar  results. 
S.  Bonder  has  labeled  this  approach  a  'fitted-parameter  analytical  model'  [Refs.  2,3:  pp. 
45,73-88]. 

The  second  method  for  developing  the  attrition-rate  coefficients  has  been 
labeled  an  'independant  analytical  model'  by  S.  Bonder  [Refs.  2,3:  pp.  47,73-88].  The 
concept  for  developing  the  coefficients  is  to  consider  a  single  firer  shooting  at  a  single 
target  (which  is  stationary  and  does  not  fire  back).  These  coefficients  are  developed 
under  perfect  conditions  with  this  method.  There  has  been  a  large  amount  of  work 
done  in  the  development  of  acceptable  attrition-rate  coefficients.  Taylor  provides  a 
good  list  of  references  in  his  book  Force-on-Force  Attrition  Modeling  [Ref.  2]. 
5.  Modern  Conditions  of  Warfare 

In  1914  when  F.  W.  Lanchester  developed  his  original  models,  he  was 
attempting  to  model  the  combat  attrition  process  so  that  he  could  learn  more  about 
the  dynamics  of  combat  under  'modern  conditions'.  The  phrase  modern  conditions 
was  an  important  one  then  as  it  is  today.  In  earlier  sections  Lanchesters  models  were 
introduced  and  the  assumptions  that  applied  to  them  were  given.  The  assumption  sets 
are  what  tailor  the  model  to  the  conditions.  Some  of  these  original  assumptions  no 
longer  seem  appropriate  to  model  what  we  would  term  'combat  under  modern 
conditions'.  The  next  section  discusses  some  of  the  shortcomings  of  the  original 
Lanchester  Models,  but  before  that  it  is  important  to  lay  out  a  solid  definition  for 
combat  under  modern  conditions. 

Today  the  conditions  on  the  combat  battlefield  are  very  complex.  First,  the 
air  and  land  battles  are  thought  of  as  one  battle  and  given  the  title  'AIR-LAND 
BATTLE'.  The  ground  forces  commander  from  the  battalion  level  on  up,  has  a 
combination  of  different  combat  units  under  his  direct  command.  An  Army  Division 
Commander  has  Infantry,  Armor,  and  Artillery  Battalions  all  assigned  to  his  command. 
The  battalions  which  are  assigned  are  not  pure'  units,  they  too  have  a  combination  of 
combat  units  assigned.  This  concept  has  been  adopted  by  the  Army  as  the  way  to 
fight  the  next  war.    The  Army  has  organized  under  the  'combined  arms'  concept  and 
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believe  that  placing  a  mix  of  supporting  combat  forces  together  on  the  battlefield  will 
be  a  force  multiplier. 

The  Navy  has  developed  much  the  same  type  of  conditions  at  sea.  The 
Carrier  Battle  Groups  (CBG)  are  combinations  of  many  types  of  naval  combat  power 
and  therefore  brings  a  variety  of  forces  together  to  fight  a  battle  at  one  point  under  the 
Ofiicer-in-Tactical-Command  (OTC). 

6.  Shortcomings  With  The  Original  Lanchester  Models 

Lanchester's  original  models  provided  a  very  simplified  model  of  modern 
warfare  which  has  served  as  the  basis  for  a  great  deal  of  further  work,  in  this  area. 
There  are  many  shortcomings  of  the  original  models  which  have  been  identified  and 
have  stimulated  much  of  the  work  done  in  the  years  following  their  development. 
Lanchesters  models  have  been  popular  because  of  their  simplicity  and  the  ease  of 
tracing  through  the  mathematical  computations.  In  this  section  I  will  discuss  some  of 
the  shortcoming  which  have  been  identified  in  applying  Lanchesters  original  laws  to 
model  today's  combat. 

First,  not  all  of  today's  combat  scenarios  can  be  neatly  placed  into  the 
'square'  or  'linear'  attrition  models.  The  most  often  referenced  and  easily  explained 
example  is  a  guerilla  ambush  scenario.  In  this  case  it  is  assumed  that  the  guerilla  force 
is  heavily  camouflaged  and  the  force  being  attacked  is  in  the  open.  It  is  apparent  that 
the  attacked  force  would  suffer  a  'squared'  attrition  rate  from  the  aimed  fire  of  the 
guerilla  force  which  has  clear  line  of  sight  of  the  entire  force  under  attack.  On  the 
other  hand,  the  defending  force  would  not  have  clear  line  of  sight  of  the  attacking 
guerilla  force  and  therefore  would  not  be  properly  modeled  by  a  'square'  attrition 
model,  but  rather  the  'linear'  attrition  model.  However,  even  this  does  not  properly 
model  the  scenario.  With  time  the  defending  force  will  be  able  to  locate  the  attacking 
force  well  enough  that  their  fire  should  be  considered  as  'aimed  fire'  and  modeled  by  a 
'squared'  attrition  model.  This  illustrates  the  dynamic  conditions  on  the  battlefield 
which  must  be  considered.   [Ref.  4] 

One  of  the  main  assumptions  of  Lanchester's  laws  is  that  the  two  opposing 
forces  are  comprised  of  homogeneous  units  (i.e.,  Infantry  vs..  Infantry).  This  is  rarely 
the  case  on  the  battlefield  today.  Today  the  doctrine  for  a  land  battle  is  the  Air-Land 
Battle  Doctrine,  and  is  written  for  combined  arms  units.  The  combined  arms  concept 
was  introduced  earlier,  but  warrants  reemphasis.  The  force  that  must  be  modeled 
today  combines  several  different  combat  units  under  one  commander  to  fight  a  battle. 


19 


An  additional  consideration  is  the  close  relationship  between  the  ground  force  and  the 
close  air  support  provided  by  the  Air  Force  which  adds  another  unit  type  into  the 
model.  The  same  type  of  unit  composition  is  currently  found  in  the  Navy  Carrier 
Battle  Groups.  The  need  for  a  model  which  allows  for  several  different  types  of  units 
(the  heterogeneous  force  model)  is  obvious.  The  use  of  heterogeneous  forces  requires 
the  model  to  consider  allocating  fire  power  on  a  variety  of  targets.  The  factors  that 
enter  in  the  allocation  process  are;  the  attrition  rate  of  each  type  of  weapon  for  each 
of  the  opposing  forces  weapons  types,  the  number  and  type  of  enemy  targets,  and  the 
relative  significance  of  each  of  the  enemy  targets. 

Another  assumption  of  the  original  Lanchester  theory  is  that  all  of  the  forces 
of  both  sides  are  committed  to  battle  at  the  beginning  and  there  are  no  reinforcements 
available.  Furthermore,  there  is  no  allowance  for  the  possibility  of  withdrawing  some 
or  all  of  a  unit's  forces.  These  factors  are  definitly  consider  by  the  commander  on  the 
battlefield  today. 

Separation  of  forces  on  the  battlefield  is  ignored  in  the  original  models.  They 
assume  that  the  units  remain  fixed  in  location  and  are  always  within  range  of  the 
opposing  forces  weapons.  To  better  model  the  combat  process,  a  movement  factor 
which  would  affect  the  weapons'  ranges  and  the  accuracy  of  each  weapon  on  the 
battlefield  would  be  required. 

The  original  model  is  designed  to  model  a  Tight  to  the  finish'  type  of  combat. 
This  definition  for  battle  termination,  of  fighting  to  the  last  man,  is  not  realistic  on 
today's  battlefield.  There  is  some  point  in  a  battle  where  the  defender  will  decide  to 
withdraw  or  surrender  before  the  entire  force  is  destroyed.   [Ref.  4] 

With  this  introduction  to  Lanchester  Equation  type  models  as  a  basis,  the 
next  chapter  will  present  the  two  models  used  for  this  research. 
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II.  THE  MODELS 

A.  GENERAL 

Two  models  were  used  for  this  research,  one  is  a  deterministic  model  and  the 
other  is  a  stochastic  computer  simulation  model.  The  deterministic  model  used  was 
developed  at  the  Naval  Postgraduate  School  by  Professor  Paul  H.  Moose  and 
Professor  John  M.  Wozencraft.  Their  model  is  a  Generalized  Lanchester  Equation 
type  model.  It  will  be  referenced  throughout  the  remainder  of  this  thesis  as  the  M/W 
model  or  the  Moose/ Wozencraft  model. 

The  stochastic  computer  simulation  model  which  was  utilized  for  this  research  is 
a  Monte  Carlo  Simulation  written  by  Professor  Don  E.  Harrison,  Jr.  at  the  Naval 
Postgraduate  School.  It  will  be  referenced  throughout  the  remainder  of  the  thesis  as 
the  C/S  model  or  the  computer  simulation  model. 

In  the  following  sections,  each  of  these  models  will  be  presented  to  provide  a 
general  understanding  of  there  characteristics. 

B.  MOOSE/WOZENCRAFT  MODEL 

Paul  H.  Moose  is  a  professor  in  the  Electrical  Engineering  Department  at  the 
Naval  Postgraduate  School  and  John  M.  Wozencraft  was  a  professor  with  the 
Electrical  Engineering  Department  prior  to  his  retirement.  Both  have  been  associated 
with  the  academic  group  for  the  Joint  Command,  Control,  and  Communications 
Curriculum,  and  have  been  interested  in  the  decision  and  control  problems  a  modern 
day  military  commander  faces.  They  wanted  to  develop  a  better  understanding  of  how 
command  and  control  decisions  affect  battle  outcome.  To  do  this,  an  analytical  model 
of  the  attrition  process  that  would  provide  an  adequate  representation  of  the  modern 
conditions  of  combat  was  needed.  Such  a  model  required  modeling  a  variety  of  forces 
as  an  aggregate  force  fighting  on  the  battlefield.   [Ref.  5] 

As  presented  in  chapter  one,  the  military  commander  in  today's  environment  has 
several  types  of  units  available  to  him  to  fight  a  battle.  The  question  of  optimum 
resource  allocation  involves  the  use  of  a  variety  of  these  units.  For  the  model  to  be 
useful  it  had  to  be  easy  to  interpret  and  understand.  Lanchester  Equation  type  models 
were  easy  to  interpret,  but  were  limited  in  their  ability  to  model  a  variety  of  forces. 
This  motivated  their  research  into  the  dynamics  of  a  Generalized  System  of  Lanchester 
Equations  of  the  type  shown  here  in  Equation  Set  2.1: 
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xi(t)=  F^x^y.t):  where  i=  1,2, N, 

(eqn2.1) 
yj(t)  =  Gj(x,yj,t):  where  j=  1,2, ,M, 

where  N  components  of  a  non-homogeneous  X-force  engage  M  components  of  a  Y- 
force.  As  shown  in  Equation  Set  2.1,  the  number  of  combatants  of  type  x-  is  a 
function  of  their  own  forces,  the  total  Y  force  (the  sum  of  the  y:  forces)  available  to 
fight,  and  time.  Similarly,  the  rate  of  change  of  a  y:  force  is  a  function  of  it's  own 
force  level,  the  total  X  force  (the  sum  of  all  x-  forces)  available  and  time. 

The  following  set  of  equations  represent  the  Moose/ Wozencraft  model: 

xi(t)  =  -uixi(t)-XaijXi(t)y-(t)-^bip-(t)  +  r{    where  i=  1,2 N, 

J                      J  (eqn  2.2) 

y-(t)  =  -vjyj(t)-XcijXi(t)yj(t)-j;dijXi(t)  +  s-      where  j  =  1,2, M, 

where 

1)  The  x-  and  y:  represent  Blue  and  Orange  forces  of  different  types 

2)  The  U|  and  V:  are  self  attrition  coefficients 

3)  The  a-:  and  c-  are  area  fire  coefficients 

4)  The  b-  and  d-  are  aimed  fire  coefficients 

5)  The  r-  and  s-  are  resupply  coefficients 

6)  x-,v:  ^  0  for  all  i  and  j 

7)  £xi=X,£y:  =  Y. 

J 

The  M/W  model  is  designed  as  an  NxM,  heterogeneous  force,  combat  attrition 
model.  The  significance  of  this  fact  is  that  it  allows  for  closer  modeling  of  the  modern 
combat  force  discussed  in  Chapter  One. 

In  the  M/W  model,  Equation  Set  2.2,  several  things  are  different  from  the 
original  Lanchester  Equation  models  presented  in  Chapter  One.  One  major  difference 
is  that  the  M/W  model  allows  for  modeling  N  type  X  forces  and  M  type  Y  forces.  The 
Moose/ Wozencraft  model  treats  the  attrition  of  a  given  force  as  one  process  where 
Lanchester's  original  work  broke  the  attrition  process  into  two  models,  aimed  and  area 
fire  models.  With  that  fact  in  mind  then,  the  rate  of  change  of  a  force  x^  is  a  function 
of  four  separate  factors;    self-attrition,  area  fire  attrition,  aimed  fire  attrition  and 
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resupply.  In  the  case  of  determining  the  attrition  to  force  x-  ,  the  self-attrition  is 
modeled  using  a  self-attrition  coefficient  -Uj  which  is  multiplied  by  the  current  x-  force 
level.  This  yields  the  attrition  to  force  Xj  which  is  attributed  to  'self  losses'  such  as 
disease  and  defections.  The  area  fire  attrition  is  modeled  similarly  to  Lanchesters 
original  model,  with  the  added  capability  to  account  for  all  the  possible  combinations 
ofy:  forces  which  can  cause  casualties  to  x-  by  means  of  area  fire.  Thus  an  attrition- 
rate  coefficient  for  each  possible  combination  of  forces  engaging  in  battle  must  be 
available.  Each  y:  force  has  an  area  fire  coefficients  -aj:  associated  with  it.  The  sum  of 
the  attrition-rate  coefficients  (a-:)  multiplied  by  the  y:  force  level  can  then  be  thought 
of  as  an  aggregate  area  fire  attrition  rate  which  multiplies  the  x-  force  level.  The  aimed 
fire  is  very  similar.  It  is  also  similar  to  the  original  Lanchester  model  with  the  added 
capabilities  to  model  NxM  force  battles.  The  aimed  fire  is  modeled  by  multipling  each 
y:  force  by  an  attrition-rate  coefficient  -bj:  and  summing  them  up.  This  is  then  the 
aggregate  aimed  fire  attrition  rate.  Finally  the  last  factor  modeled  is  the  resupply  of 
each  x-  force  with  a  resupply  variable  r{.  The  sum  of  the  self-losses,  area  fire  losses, 
aimed  fire  losses  and  resupply  is  equal  to  the  change  in  the  x-  force. 

The  model  determines  the  change  in  a  y:  force  in  exactly  the  same  manner.  The 
attrition  rate  coefficients  are  different  of  course,  but  the  process  is  the  same.  The  self- 
attrition  coefficient  for  y:  is  -V:.  The  area-fire  attrition  coefficients  are  given  as  -Cj:  and 
the  aimed  fire  coefficients  are  given  as  -djj.   The  resupply  variable  for  y:  is  S:. 

The  relative  usefulness  of  this  model  for  modeling  the  combat  attrition  process 
on  today's  battlefield  is  apparent  from  this  discussion.  First,  it  provides  the  capability 
to  model  a  battle  involving  two  heterogeneous  forces  which  we  saw  is  necessary  to 
analyze  the  optimum  force  allocation  issue.  Secondly,  it  is  a  deterministic  model,  so  it 
is  easy  to  understand  and  interpret.  It  has  many  similarities  to  the  original  Lanchester 
Equation  models  which  are  widely  studied  and  understood. 

The  results  of  Professor  Moose  and  Professor  Wozencraft's  research  into  the 
dynamical  properties  of  this  model  are  provided  in  a  paper  which  they  are  submitting 
to  the  Military  Operations  Research  Society  for  publication  titled  Characteristic 
Trajectories  of  Generalized  Lanchester  Equations  [Ref.  6]. 

C.       STOCHASTIC  COMPUTER  SIMULATION  MODEL 

The  stochastic  model  is  a  Monte  Carlo  Simulation  designed  to  test  the 
Moose/ Wozencraft  model.    The  basis  of  this  model  was  designed  by  Don  E.  Harrison, 
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Jr.,  a  professor  at  the  Naval  Postgraduate  School.  Professor  Harrison  was  asked  to 
design  a  stochastic  model  of  the  attrition  process  which  included  the  same  factors  as 
the  Moose/ Wozencraft  model.  The  combination  of  his  efforts  and  some  additional 
code  to  run  the  program  for  a  given  number  of  replications  and  calculate  statistical 
data  resulted  in  the  model  shown  in  Appendix  B  entitled  Computer  Simulation  model 

(C/S). 

The  C/S  model  will  be  presented  in  the  following  sections: 

1)  Input  parameters 

2)  The  Combat  Cycle 

3)  The  Output 

The  material  in  each  section  is  presented  in  a  condensed  users  manual  form  to  provide 
the  reader  with  an  understanding  of  the  model  and  the  capability  to  exercise  the  model 
if  desired. 

1.  Input  parameters 

The  computer  code  is  written  in  FORTRAN77  and  uses  an  exec  file  to  define 
all  the  input  files,  output  files,  load  and  execute  the  program. 

The  input  file  contains  all  the  required  input  data.  Figure  2.1  shows  the 
contents  and  format  of  the  input  data  file.  The  first  line  of  the  file  is  a  heading  which 
is  read  by  the  program  and  allows  the  user  to  label  and  keep  track  of  which  test  data 
set  is  being  utilized.  Next  the  data  required  for  the  x-  forces  is  given.  On  the  second 
line  the  number  of  types  of  X  forces  is  given  (the  variable  NTYPX  in  the  program), 
this  corresponds  to  the  N  value  in  the  M/W  model  which  specifies  the  number  of  types 
of  X  forces.  Then  the  initial  force  level  for  each  x-  force  is  given,  their  locations,  their 
corresponding  self-attrition  coefficients  (u^)  and  the  resupply  variables  (r-).  Next,  the 
attrition-rate  coefficients  are  provided  in  matrix  form.  The  aimed  fire  attrition-rate 
coefficients  'd-'  are  given  in  lines  three  through  five.  Area  fire  attrition- rate 
coefficients  'c-:'  are  given  in  lines  six  through  eight.  Then  the  required  information  for 
the  y:  forces  is  listed.  Once  again,  the  number  of  type  of  y:  forces  is  given  (the 
NTYPY  variable  in  the  program)  which  corresponds  to  M  in  the  M/W  model.  Then 
each  y:  initial  force  level,  location,  self  attrition-rate  coefficients  'v:'  and  resupply 
variables  's:'  are  given.  Lines  ten  through  twelve  contain  the  aimed  fire  coefficients 
'b-:'  and  lines  thirteen  through  fifteen  contain  the  area  fire  coefficients  'a-'.  Line 
sixteen  of  the  file  contains  seed  numbers  for  the  pseudo-random  number  generators. 
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1  01/16/87  TEST  OF  (2X2),  1  UNSTABLE  ROOT  S 

2  NSTEP 

3  NTYPX  NTX1  NTX2  NTX3  loci  loc2  loc3  Uj  u2  u3  rj  r2  r3 

4  dnd12d13 

5  d21d22d23 

6  d31d32d33 

7  cll  c12c13 

8  c21  c22  c23 

9  c31  c32c33 

10  NTYPY  NTY1  NTY2  NTY3  loci  loc2  loc3  \{  v2  v3  sj  s2  s3 

11  bn  b12bn 

12  b21  b22  b23 

13  b31  b32b33 

14  au  a12a13 

15  a21  a22  a23 

16  a31  a32  a33 

17  DSEED  DSEED  DSEED 


Figure  2.1     Input  Data  File. 

The  program  is  designed  to  be  interactive  before  it  begins  the  combat  cycle 
which  allows  the  user  to  change  the  heading  of  the  input  file  if  so  desired.    This  also 
enables  the  user  to  ensure  that  the  input  file  is  the  one  desired. 
2.  The  Combat  Cycle 

The  combat  cycle  is  a  stochastic  process.  During  the  combat  cycle  the 
combatants  are  chosen  at  random  and  the  probability  of  self-loss,  aimed  fire  loss,  and 
area  fire  losses  are  all  tested  by  use  of  a  pseudo-random  number  between  0.0  and  1.0. 
The  cycle  begins  by  randomly  selecting  which  force,  X  or  Y,  will  fire.  Once  the  force 
has  been  selected  then  a  combatant  from  the  force  is  selected.  For  the  purpose  of 
presentation,  say  that  Xj(l)  (that  is  combatant  one  of  force  type  Xj)  is  selected  to  be 
the  firer,  then  the  first  thing  that  occurs  in  the  cycle  is  a  check,  to  see  if  combatant 
Xj(l)  is  a  self  loss.   This  is  accomplished  by  use  of  a  random  number  (between  0.0  and 
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1.0)  which  is  compared  against  the  appropriate  self  attrition  coefficient  Ui.  If  the 
combatant  is  not  a  self  loss  then  a  Y  force  combatant  is  selected  at  random  as  a  target, 
say  y7(2)  (the  second  combatant  of  force  type  y^\  as  an  example.  Once  the  firer  and 
target  have  been  identified,  the  next  step  tests  the  effectiveness  of  the  'shot'  by  using  a 
random  number  to  test  whether  the  target  is  killed.  This  process  involves  first 
determining  if  the  target  is  killed,  which  is  done  by  testing  whether  a  random  number  is 
less  than  the  total  kill  probablity  (which  is  the  sum  of  the  aimed  fire  and  area  fire 
attrition-rate  coefficients,  or  their  kill  probabilities).  If  the  random  number  is  less,  then 
the  target  has  been  killed  and  if  not  the  shot  missed.  If  the  target  was  killed  then  this 
same  process  is  followed  to  determine  whether  it  was  a  loss  due  to  aimed  fire  or  area 
fire.  Required  data  is  tabulated  on  the  number  of  combatants  of  type  y:  which  are 
killed  by  type  x-  forces  and  each  force  is  resupplied.  The  combat  cycle  is  repeated  in 
this  manner  until  all  combatants  have  had  an  opportunity  to  fire  during  each  timestep. 
This  process  is  continued  for  each  timestep  until  one  of  two  conditions  occurs: 

1)  One  force  reaches  a  break  point  which  has  been  defined  as  NXSTOP  and 
NYSTOP.  This  simulates  the  level  of  attrition  a  commander  will  suffer  before 
pulling  back  or  surrendering  rather  than  fighting  to  the  last  man. 

2)  The  combat  cycle  has  been  repeated  for  a  given  number  of  timesteps,  which  is 
specified  in  the  input  data  file. 

3.  The  Output 

The  C/S  model  was  designed  to  provide  several  output  files  for  data  analysis. 
The  model  generates  two  files  which  provide  a  very  detailed  timestep-by-timestep 
recording  of  the  combat  cycle  outcome.  These  files  require  a  great  deal  of  storage 
space  and  therefore  during  the  interactive  portion  of  the  program  the  user  is  asked  if 
they  want  these  files  printed  in  their  complete  form.  It  is  not  recommended  that  they 
be  printed  if  the  user  desires  multiple  replications  because  the  storage  space  required  is 
too  large.  One  of  these  files  provides  a  detailed  listing  of  each  timestep  results 
including  a  summary  of  the  force  levels,  the  number  of  losses  by  each  type  of  fire  and 
the  number  of  combatants  that  were  resupplied.  An  example  of  this  file's  output  is 
shown  in  Figure  2.2.  The  other  file  contains  similar  listings  with  the  number  of  dead 
combatants  of  each  type  provided  as  well.  An  example  of  this  file's  output  is  shown  in 
Figure  2.3. 

The  remainder  of  the  files  are  designed  to  provide  the  force  levels  of  each  force 
(i.e.,  Xj  and  yj)  for  each  timestep.    This  provides  the  data  which  can  be  used  to  plot 
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02/02/87  TEST  OF  (3X3), 

ft  X  TYPES  NUMBER  EACH    LOCATIONS      SELF-LOSS       RE -SUPPLY 

3   200   150   175   1.0   1.0   1.0   0.005   0.010   0.020    4    4    11 

KILL  PROBABILITIES  FOR  FORCE  X(  I )  SHOOTS  Y(J) 

PK(1,1)  PK(1,2)  PK(1,3)  PK(2,1)  PK<2,2)  PK(2,3)  PK(3,1)  PK<3,2)  PK(3,3) 
0.06000  0.02000  0.02000  0.00000  0.00000  0.00000  0.03000  0.01000  0.02000 


AREA  FIRE  COEFFICIENTS  FORCE  X(I)  SHOOTS  Y(J) 

CC(1,1)  CC(1,2)  CC(1,3)  CC(2,1)  CC(2,2)  CC(2>3)  CC(3,1)  CC(3,2)  CC<3,3) 
0.00000  0.00000  0.00000  0.00000  0.00030  0.00020  0.00000  0.00000  0.00010 

ft  Y  TYPES  NUMBER  EACH    LOCATIONS      SELF-LOSS       RE-SUPPLY 

3   150   125   225   1.0   1.0   1.0   0.003   0.020   0.030    4    5   14 

KILL  PROBABILITIES  FOR  FORCE  Y(I)  SHOOTS  X(J) 

PK<1,1)  PK(1,2)  PK<1,3)  PK(2,1)  PK(2,2)  PK(2,3)  PK(3,1)  PK(3,2)  PK(3,3) 
0.05000  0.01000  0.05000  0.00000  0.00000  0.00000  0.02000  0.01000  0.03000 


AREA  FIRE  COEFFICIENTS  FORCE  Y(I)  SHOOTS  X(J) 

AA(1,1)  AA(1,2)  AA(1,3)  AA(2,D  AA(2,2)  AA(2,3l  AA(3,1)  AA(3,2)  AA(3,3) 
0.00000  0.00000  0.00000  0.00000  0.00020  0.00030  0.00000  0.00000  0.00010 

DOUBLE  PRECISION  RANDOM  VARIANT  SEEDS  INTEGER  SEED  VALUES 

0. 345215790+09  Q.45635761D+08  0.89342761D+08  1209308051  1194034571  1196766258 

ft   NTXS  NTX  NX1   NX2   NX3  RX1  RX2  RX3  NTYS  NTY  NY1   NY2   NY3  RY1  RY2  RY3 
+   0   525  525   200   150   175   4   4   11   500   500   150   125   225   4   5   14 


TIMESTEP    1 

AIMED  FIRE  KILLS  OF  Y"S  BY  X  FORCE 

D(l,l)   D(l,2)   D(l,3)   0(2,1)   0(2,2)   D(2,3)   0(3,  II   0(3,2)   D(3,3) 

302000110 

AIMED  FIRE  KILLS  OF  X"S  BY  Y  FORCE 

B(l,l)   B(l,2)   B(l,3)   B(2,l)   B(2,2)   B(2,3)   B(3,l)   B(3,2)   B(3,3) 

211000202 

AREA   FIRE   KILLS  OF  Y"S     BY  X   FORCE 

C(l,l)     C(l,2)     C(l,3)     C(2,l)     C(2,2)      C(2,3)      C(3,l)     C(3,2)      C(3,3) 

000057002 

AREA   FIRE   KILLS  OF  X"S     BY  Y   FORCE 

A(l,l)      All, 2)      All, 3)      A(2,l)      A(2,2)      A(2,3)      A(3,l)      A(3,2)      A13.3) 

000025003 

ft  NTXS  NTX  NX1  NX2  NX3  DX1  0X2  0X3  NTYS  NTY   NY1  NY2   NY3  DY1  DY2  DY3 
+    1   525  525   200   150   175   5   5   13   500   500   150   125   225   4    9   20 


Figure  2.2    Example  File  6  Output. 
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TOTAL  X  FORCE  LOSS  WAS   23,  RESUPPLY  WAS:    4   4   11 
AIMED  FIRE  LOSSES:    4    13 
AREA  FIRE  LOSSES:    0    2   8 

SELF  LOSSES:    12    2 
TOTAL  Y  FORCE  LOSS  WAS   33,  RESUPPLY  WAS:    4   5   14 
AIMED  FIRE  LOSSES:    4    12 
AREA  FIRE  LOSSES:    0   5    9 

SELF  LOSSES:    0    3    9 
1   525  521   199   149   173   4   4   11   500   490   150   121   219   4   5   14 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 

AA(1,1)  AA(1,2)  AA(1,31  AA(2,1)  AA(2,2)  AA(2,3)  AA(3,1)  AA(3,2)  AA(3,3) 
0.00000  0.00000  0.00000  0.00000  0.00011  0.00023  0.00000  0.00000  0.00008 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 

BB(1,1)  BB(1,2)  BB(1,3)  BB(2,1)  BB(2,2)  BB(2,3)  BB(3,1)  BB(3,2)  BB(3,3) 
0.03500  0.02333  0.02000  0.00000  0.00000  0.00000  0.02333  0.00000  0.02667 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

CC(1,1)  CC<1, 2)  CC(1,3)  CC(2,1)  CC(2,2)  CC(2,3)  CC(3,1)  CC(3,2)  CC(3,3) 
0.00000  0.00000  0.00000  0.00000  0.00027  0.00021  0.00000  0.00000  0.00005 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

DDd.li    DD(1,2)    DD(1,3)   DD(2,1)   DD(2,2)   DD(2,3)    DD(3,1)   DD(3,2)   DD(3,3) 
0.05000   0.00000   0.02222   0.00000   0.00000   0.00000   0.01905   0.02286   0.00000 

STOPPED  BY  PROGRAM  AT  NSTOP  =     1 


Figure  2.2    Example  File  6  Output,  (cont'd.) 

force  level  trajectories,  which  allows  for  ease  in  performing  a  comparison  of  the  two 
model's  results  for  a  given  input  set. 

Now  that  each  of  the  models  has  been  presented,  chapter  three  will  discuss  the 
design  of  the  experiment,  and  provide  samples  of  the  output  files  used  for  plotting  the 
force  level  trajectories. 
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02/02/87  TEST  OF  (3X3), 
«  NTXS  NTX  NX1   NX2   NX3  DX1  DX2  0X3  NTYS  NTY   NY1   NY2   NY3  DY1  DY2  DY3 
+   0   525  525   200   150   175   0    0   0   500   500   150   125   225   0   0   0 


D,B,C,A  BY  ROWS  IN  DESCENDING  ORDER 

13       0 

2       0 

0 

0 

1 

1 

0 

12       1 

1       0 

0 

0 

2 

0 

2 

10       0 

0       0 

5 

7 

0 

0 

2 

10       0 

0       0 

2 

5 

0 

0 

3 

TOTAL  X  FORCE  LOSS  WAS 

23,  RESUPPLY  WAS: 

4 

4 

11 

AIMED  FIRE  LOSSES: 

4   1   3 

AREA  FIRE  LOSSES: 

0    2   8 

SELF  LOSSES:    12    2 

TOTAL  Y  FORCE  LOSS  WAS 

33,  RESUPPLY  WAS: 

4 

5 

14 

AIMED  FIRE  LOSSES: 

4   12 

AREA  FIRE  LOSSES: 

0   5    9 

SELF  LOSSES:    0   3    9 

1   525  521   199   149   173 


5   13   500  490   150   121   219 


9   20 


A,B,C,D  BY  ROWS  IN  DESCENDING  ORDER 
AV  0.00000  0.00000  0.00000  0.00000  0.00011  0.00023  0.00000  0.00000  0.00008 
AV  0.03500  0.02333  0.02000  0.00000  0.00000  0.00000  0.02333  0.00000  0.02667 
AV  0.00000  0.00000  0.00000  0.00000  0.00027  0.00021  0.00000  0.00000  0.00005 
AV  0.05000  0.00000  0.02222  0.00000  0.00000  0.00000  0.01905  0.02286  0.00000 
STOPPED  BY  PROGRAM  AT  NSTOP  =     1 


Figure  2.3     Example  File  7  Output. 
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III.  RESEARCH  METHODOLOGY 

A.  OBJECTIVE 

The  primary  objective  of  this  research  is  to  determine  whether  a  deterministic 
combat  attrition  model  is  a  good  approximation  for  the  mean  result  of  a  stochastic 
model  of  the  combat  attrition  process.  With  that  objective  in  mind,  an  experiment  was 
designed  to  preform  a  comparison  of  two  such  models.  The  models  used  in  this 
experiment  are  the  Moose/Wozencraft  model  (M/W  model),  which  is  a  deterministic 
model,  and  a  Monte  Carlo  Simulation  (C/S  model)  which  were  introduced  in  chapter 
two.  A  secondary  objective  is  to  see  if  the  two  models  behave  with  similar 
characteristics  as  the  results  of  Professor  Moose  and  Professor  Wozencraft's  research 
[Ref.  6]. 

B.  EXPERIMENTAL  DESIGN 

The  experiment  is  designed  to  run  the  C/S  model  for  ninety-nine  timesteps  over 
thirty  replications,  to  generate  force  level  data  and  calculate  the  resulting  average 
attrition  rate  coefficients.  These  average  coefficients  are  then  used  as  the  input 
parameters  to  the  M/W  model  which  then  provides  the  calculated  force  levels  per 
timestep.  This  force  level  data  is  used  for  analysis.  The  method  for  analysis  which  has 
been  chosen  is  to  plot  the  resulting  force  level  data  over  time  to  provide  for  an  easier 
comparison  than  sorting  through  the  enormous  amount  of  data  generated  for  each  test 
case. 

This  experiment  was  conducted  for  cases  involving  lxl,  1x2,  2x2,  and  3x3  force 
level  scenarios.  Examples  of  the  data  collected  and  the  resulting  force  level  trajectory 
plots  are  provided  in  the  next  section. 

1 .  Model  Verification 

Before  beginning  the  experiment  using  these  two  models,  each  of  the  models 
was  exercised  and  tested  heavily  to  verify  that  they  ran  properly  and  provided  the 
desired  output  as  they  were  designed.  A  brief  description  of  the  verification  process  for 
each  model  will  be  given. 

a.   Computer  Simulation  Model  (C/S  model) 

The  first  step  in  verifying  the  C/S  model  involved  a  procedure  of  several 
runs  using  different  input  data  sets  to  test  the  program  to  ensure  that: 
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1)  The  input  data  and  program  were  compatible  and  the  input  variables  are 
properly  read. 

2)  All  the  computational  code  was  performing  properly  (i.e.,  the  average  attrition 
rate  coefficients  and  the  average  force  levels). 

3)  That  each  of  the  test  conditions  was  functioning  as  designed  (i.e.,  the 
condition  for  ending  at  the  specified  force  levels  XSTOP  and  YSTOP). 

4)  That  the  output  was  presented  as  desired  for  the  step-by-step  analysis  and  the 
plotting  data. 

5)  Test  the  resulting  average  attrition-rate  coefficients  with  the  input  probabilities 
of  kill  which  correspond  to  these  coefficients. 

During  this  portion  of  the  research  it  was  interesting  to  note  that  the  C/S 
model  was  designed  to  allow  each  shooter  one  shot  per  time  interval,  but  the  shooter 
could  select  a  target  to  fire  at  which  had  already  been  killed  during  that  time  interval. 
This  fact  seemed  to  be  a  fairly  realistic  condition  that  could  occur  on  the  battle  field 
where  two  combatants  fire  at  the  same  target  instantaneously  or  within  a  fraction  of  a 
second  of  each  other.  Therefore,  this  condition  was  left  in  the  model. 
b.   Moose/  Wozencraft  Model  (A//  W  model) 

A  computer  program  was  designed  to  run  the  M/W  model  to  allow  for  an 
easy  comparison  of  the  two  models.  By  utilizing  a  computer  program  it  allowed  for 
the  same  timestep  by  timestep  analysis  as  the  C/S  model. 

After  designing  the  FORTRAN  program,  shown  in  Appendix  D,  to  run  the 
Moose/Wozencraft  model  (M/W  model)  it  was  tested  to  verify  that  the  following 
conditions  were  met: 

1)  The  input  data  and  the  program  were  compatible  and  all  variables  were  read 
properly. 

2)  All  computational  code  was  performing  as  designed. 

3)  The  resulting  output  files  provided  the  desired  data  for  analysis.  Since  the 
M/W  model  is  a  deterministic  model,  the  same  output  will  be  obtained  each 
time  for  the  same  set  of  input  data  parameters. 

Each  of  the  conditions  listed  above  was  verified  for  simple  input  data  sets  for  each  of 
the  cases. 

2.  Data  Generation  Procedure 

The  procedure  followed  in  conducting  this  experiment  was  a  three  step  process 
for  each  unique  case  tested.  The  first  step  involved  running  the  C/S  model  for  ninety- 
nine  timesteps  over  thirty  replications  which  would  generate  the  following  data  files: 
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1)       The  calculated  average  attrition-rate  coefficients  which  are  later  used  as  input 
to  the  M/W  model.   An  example  file  is  shown  in  Figure  3.1. 


01/27/87  TEST 

OF  (1X1) 

AVERAGE 

ATTRITION  COEFFICIENTS 

FOR  Y  FORCE 

AA(1,1) 

AA(1,2) 

AA<1,3) 

AA(2>1)  AA(2>2)  AA(2,3) 

AA(3,1)  AA(3,2) 

AA(3,3) 

0.00012 

0.00000 

0.00000 

0.00001  0.00000  0.00000 

0.00000  0.00000 

0.00000 

AVERAGE 

ATTRITION  COEFFICIENTS 

FOR  Y  FORCE 

BB(1,1) 

BB<1,2) 

BBC  1,3) 

BB(2,1)  BB(2,2)  BB<2,3) 

BB<3,1)  BB(3,2) 

BB(3,3) 

0 . 04857 

0.00000 

0.00000 

0.00142  0.00000  0.00000 

0.00000  0.00000 

0.00000 

AVERAGE 

ATTRITION  COEFFICIENTS 

FOR  X  FORCE 

CC(1,1) 

CC(1,2) 

CC(1, 3) 

CC(2,1)  CC(2,2)  CC(2,3) 

CC(3,1)  CC(3,2) 

CC(3,3) 

0.00027 

0.00000 

0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000 

0.00000 

AVERAGE 

ATTRITION  COEFFICIENTS 

FOR  X  FORCE 

DDI1>1) 

DD(1,2) 

DD(1,3) 

00(2,1)  DD(2,2)  DD(2,3) 

DD(3,1)  00(3,21 

DD(3,3) 

0.08648 

0.00000 

0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000 

0.00000 

STOPPED 

BY  PROGRAM  AT  NSTOP  =    99 

Figure  3.1     Example  Average  Attrition-rate  Coefficients  File  For  C/S  Model. 

2)  A  file  containing  the  average  force  level  for  each  force,  the  average  force  level 
plus  one  standard  deviation,  and  the  average  force  level  minus  one  standard 
deviation  is  provided  for  the  aggregate  force  levels  X  (£x:)  and  Y  (]Ty:)  per 
timestep  as  well  as  for  each  of  the  individual  units  x-  and  y:.  An  example  file 
is  shown  in  Figure  3.2. 

3)  The  average  force  level  and  standard  deviation  of  each  force  level  per  time 
step.   An  example  of  these  results  is  shown  in  Figure  3.3. 

The  second  step  involves  using  the  average  attrition-rate  coefficients  calculated 
by  the  C/S  model  as  input  parameters  to  the  M/W  model.  The  output  generated  by 
the  M/W  model  for  each  set  of  input  data  is  the  following: 

1)       The  force  level  for  each  of  the  forces  x-  and  y:  as  well  as  the  total  force  levels 

X  (£xj)  and  Y  (Yyp  for  each  time  step.    An  example  of  this  output  file  is 

shown  in  Figure  3.4. 
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01/16/87 
T    xj 

1  203.8 

2  207.1 

3  210.7 

4  212.0 

5  215.2 

TEST  OF  (2X2), 

yj       xj  +  ff 

313.9      209.2 
322.0      213.8 
327.4      219.1 
329.3      223.0 

327.8      225.7 

1  UNSTABLE  ROOT 

Vi+or      Xi-<r     yi-cr 

320.2      198.3      307.6 
331.5      200.4      312.4 
336.9      202.3      317.9 
340.1      201.0      318.4 
340.4      204.7      315.2 

97  250.1 

98  251.3 

99  252.5 

284.8      259.8 
283.7      261.9 
283.0      262.9 

296.2      240.3 
294.7      240.7 
294.9      242.1 

273.4 
272.8 
271.2 

Figure  3.2     Example  Average  XjVi  Force  Level  File  For  C/S  Model. 


01/27/87   TEST  OF  (1X2),  < 
T    x,         (Txl         yj 

1  493.43       6.64     170.37 

2  492.27       9.54     149.47 

3  496.23      11.75     132.63 

4  502.47      12.29     119.67 

5  510.30      16.21     109.67 

CASE  5 

°yl 

5.16 
6.24 
6.63 
6.35 
6.38 

95  939.63 

96  939.80 

97  940.13 

98  941.07 

99  942.87 

21.95 
25.08 
26.22 
25.40 
26.01 

67.10 
67.20 
68.20 

68.77 
68.43 

4.42 
4.43 
4.37 
4.44 
4.14 

Figure  3.3     Example  File  of  Force  Level  and  Standard  Deviation   For  C/S  Model. 

The  third  step  involves  plotting  the  force  level  data  generated  by  each  of  the 
models  to  perform  the  comparison.  The  plots  which  were  utilized  for  the  comparison 
are  listed  below  with  examples  provided: 

1)       A  plot  of  the  total  X  force  vs.  Total  Y  force  level  trajectories  for  both  models 

is  done  first.    This  plot  provides  for  a  macro  level  comparison  of  the  resulting 
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01/27/87  TEST  OF  (1x2),  CASE  4 
T     x,  y. 


0 

200.00 

300.00 

1 

196.29 

299.19 

2 

192.88 

298.75 

3 

189.73 

298.61 

4 

1S6.82 

298.76 

5 

184.12 

299.13 

95 

91.69 

472.51 

96 

90.92 

474.95 

97 

90.14 

477.41 

98 

89.37 

479.89 

99 

88.60 

482.39 

Figure  3.4    Example  K-fij  F°rce  Level  Data  File  For  M/W  Model. 

force  levels  over  time  (i.e.,  the  Force  Commanders  level).  This  plot  is  two 
dimensional,  plotting  the  X  force  level  vs..  Y  force  level.  The  time  variable  is 
shown  by  marking  the  force  level  every  fifth  timestep.  This  plot  will  be  used 
to  analyze  the  two  aggregate  force  models  results.  An  example  of  the  total 
force  level  trajectory  plot  is  shown  in  Figure  3.5. 
2)  Plots  of  the  x-  force  level  vs..  time  and  the  y:  force  levels  vs.  time  for  both  the 
C/S  and  the  M/W  models  are  given.  This  provides  a  comparison  of  the  battle 
outcome  for  the  two  models  at  the  micro  level.  An  example  of  the  x-  (or  y:) 
plot  is  shown  in  Figure  3.6. 
The  plots  introduced  in  this  chapter  are  used  extensively  for  conducting  the 

comparison  of  the  two  models.   Chapter  Four  will  present  the  results  of  the  experiment 

and  discuss  the  comparisons. 
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Figure  3.5    Example  Total  Force  Level  Trajectory  Plot. 
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XI  FORCE  LEVEL  TRAJECTORY 
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Figure  3.6    Example  plot  of  x^  force  level  trajectories. 
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IV.  EXPERIMENTAL  RESULTS 

A.  INTRODUCTION 

The  primary  objective  of  this  experiment  was  to  test  whether  the  deterministic 
Moose/ Wozencraft  model  was  a  good  approximation  (or  prediction)  for  the  mean 
result  of  the  stochastic  computer  simulation.  The  secondary  objective  for  this 
experiment  was  to  test  whether  the  two  models'  force  level  trajectories  behaved  in  the 
expected  manner  described  by  Professors  Moose  and  Wozencraft.  The  results  of  ten 
cases  are  presented  graphically  showing  the  two  models'  results  with  a  bracket  of  plus 
and  minus  one  standard  deviation  as  a  measure  of  uncertainty.  Each  case  will  be 
discussed  in  the  context  of  the  primary  and  secondary  objectives  in  this  chapter  and 
their  results  shown  graphically  in  Appendix  E. 

B.  THE  FOUR  MODES  OF  BEHAVIOR 

Before  proceeding  with  the  presentation  of  the  results,  it  is  important  to  briefly 
explain  the  four  trajectory  behavior  modes  which  will  be  referenced  in  each  of  the 
cases.  Professors  Moose  and  Wozencraft  found  that  the  force  level  trajectories 
displayed  one  of  four  distinct  modes  of  behavior.  A  mode  of  behavior  was  the  result  of 
the  existence,  or  nonexistence,  of  equilibrium  points  in  the  state  space  (positive  first 
quadrant).  The  four  modes  related  to  the  existence  of  one  stable  equilibrium,  one 
unstable  equilibrium,  two  equilibria  (where  one  is  stable  and  one  is  unstable),  or  no 
equilibria.   These  four  modes  of  behavior  are  illustrated  in  Figure  4.1.   [Ref.  6] 

The  expected  behavior  for  any  force  level  trajectory  when  one  stable  equilibrium 
point  exists  in  the  state  space  is  shown  in  Figure  4.1  (a).  The  trajectory  will  be 
attracted  to  the  stable  equilibrium  and  remain  there,  since  dx/dt=dy/dt  =  0  at  the 
stable  equilibrium. 

The  expected  behavior  for  any  force  level  trajectory  when  one  unstable 
equilibrium  point  exists  in  the  state  space  is  shown  in  Figure  4.1  (b).  As  the 
trajectories  near  the  equilibrium  they  are  repelled  or  turned  away  from  the  equilibrium 
point.  The  presence  of  an  unstable  equilibrium  creates  a  division  of  the  state  space 
into  two  regions.  The  dividing  line  was  called  a  separatrix.  In  scenarios  which  exhibit 
this  mode  of  behavior  the  outcome  depends  on  which  side  of  the  separatrix  the  initial 
force  levels  begin. 
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Figure  4.1     The  Four  Trajectory  Modes  Of  Behavior. 
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The  expected  behavior  for  any  force  level  trajectory  when  two  equilibria  exist 
(one  stable  and  one  unstable)  is  shown  in  Figure  4.1  (c).  Once  again  the  unstable 
equilibrium  creates  a  separatrix  dividing  the  state  space  into  regions.  The  stable 
equilibrium  will  attract  trajectories  as  discussed  earlier.  This  particular  mode  of 
behavior  has  two  possible  outcomes  for  a  battle,  one  of  the  forces  can  win  or  they  will 
fight  to  a  draw.   The  outcome  depends  on  the  initial  force  levels  for  each  force. 

The  expected  behavior  for  any  force  level  trajectory  when  no  equilibrium  point 
exists  is  shown  in  Figure  4.1  (d).  The  trajectories  will  approach,  and  eventually  cross 
one  of  the  axes,  which  results  in  a  force  being  eliminated. 

C.       DISCUSSION 

This  numerical  experiment  was  designed  to  test  a  combat  scenario  involving 
different  opposing-force  composition.  (N„  x  My)  cases  involving  Nx  types  of  X  forces 
vs.  M  types  of  Y  forces  are  tested  in  this  experiment.  An  example  is  a  1x2  case  which 
involves  one  X  type  force  vs.  two  Y  type  forces.  The  input  data  sets  were  chosen  to 
test  the  two  models  for  the  four  distinct  modes  of  behavior  that  Professors  Moose  and 
Wozencraft  identified  [Ref.  6].  The  computed  results  are  presented  in  sections  which 
are  dedicated  to  a  specific  (X  x  M)  scenario. 

Appendix  E  contains  the  input  data  sets  used  for  each  model,  and  the  resulting 
force  level  trajectory  plots  for  each  of  the  cases.  In  each  case  the  input  data  for  the 
C/S  model  was  developed  to  test  one  of  the  modes  of  behavior  and  the  M/W  model 
input  data  were  the  calculated  attrition  rate  coefficients  calculated  by  the  C/S  model. 
Several  trajectory  plots  were  done  for  each  of  the  cases.  There  is  a  plot  of  the  total 
aggregate  force  level  trajectory  which  compares  the  total  X  force  vs.  the  total  Y  force. 
This  aggregate  trajectory  shows  the  battle  outcome.  An  example  of  this  plot  is  shown 
for  case  six  in  Figure  4.2. 

The  result  of  the  computer  simulation  model,  which  is  the  mean  of  thirty 
replications,  is  always  shown  as  a  solid  line  with  squares  marking  the  force  level  every 
fifth  time  step.  The  Moose/Wozencraft  model  results  are  shown  as  a  dashed  line  with 
circles  marking  the  force  level  every  fifth  time  step.  Additionally,  two  solid  lines  are 
shown  with  no  markings.  These  are  the  curves  which  show  plus  and  minus  one 
standard  deviation  from  the  C/S  model  results.  They  serve  as  a  measure  of  uncertainty 
for  the  comparison.  The  aggregate  force  level  plot  provides  a  comparison  of  the  two 
models  total  force  level  results  to  determine  whether  they  do  have  the  same  general 
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Figure  4.2    Total  Force  Level  Trajectory  For  Case  Eight. 
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behavior  and  whether  the  deterministic  model  is  a  good  approximation  for  the  result  of 
the  stochastic  model.  Each  of  the  aggregate  forces  (i.e.,  X  =  ]Tx-)  are  plotted  as  a 
function  of  time  to  provide  a  clearer  comparison  of  the  two  models  results  as  the  battle 
progresses.  An  example  of  these  plots  are  shown  in  Figures  4.3  and  4.4.  This  provides 
another  comparison  between  the  two  models'  results  for  each  aggregate  force.  In 
addition  to  the  aggregate  force  level  plots,  the  results  for  each  of  the  separate  units  are 
plotted  as  a  function  of  time.  Examples  of  these  plots  are  provided  in  Figures  4.5 
through  4.8.  This  provides  a  comparison  of  the  two  models  at  the  individual  unit  level. 
The  method  of  comparison  used  for  this  analysis  was  a  graphical  presentation  of 
the  force  level  trajectories  rather  than  an  extensive  numerical  analysis  of  the  resulting 
data.  Each  case  will  be  presented  in  a  similar  manner.  The  cases  will  be  briefly 
described  and  then  a  list  of  observations  given.  Specific  comments  are  made  about 
whether  the  M/W  trajectory  was  close  to  the  C/S  trajectory  and  whether  the  expected 
behavior  was  observed. 

D.       (IX  1)  CASE  COMPARISONS 

The  experiment  was  conducted  for  cases  one  through  four  using  one  X  force  and 
one  Y  force,  which  are  labeled  as  lxl  cases.  The  results  of  each  of  these  cases  are 
shown  in  Figures  E.l  through  E.12  in  Appendix  E.  The  force  level  trajectories  for  the 
M/W  model  were  very  close  to  the  mean  force  level  trajectories  for  the  C/S  model  in  all 
four  cases. 

1.  Case  One 

Case  one  was  designed  to  test  the  two  models  involving  a  lxl  scenario  for  the 
mode  of  behavior  where  no  equilibrium  exists.  The  following  observations  are  noted 
from  the  results  shown  in  Figures  E.l  through  E.3: 

a)  The  resulting  force  level  trajectories  for  the  M/W  model  and  the  C/S  model 
were  very  similar,  both  displaying  the  same  trajectory  shape. 

b)  The  M/W  model  trajectory  was  well  within  the  one  standard  deviation 
boundaries,  for  the  complete  battle.  In  fact,  the  two  curves  over  lapped 
during  portions  of  the  battle.  This  case  resulted  in  the  X  force  winning  the 
battle. 

2.  Case  Two 

Case  two  was  designed  to  test  the  mode  when  one  stable  equilibrium  exists  in 
the  vicinity  of  (200,200).  The  following  observations  are  noted  from  the  results  shown 
in  Figures  E.4  through  E.6: 
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Figure  4.3    X  Force  Level  Trajectory  Over  Time  For  Case  Eight. 


42 


Y  FORCE  LEVEL  TRAJECTORY 


o 

.  o 

-J  «n 

W 
-J  o 

o 

i-u  o 
o 


o 
«o  - 
n 


o 

o 


o 


W 


LEGEND 
o    C/S  MODEL 


MSDEE 

c7s'PLU$  1  SD" 

7S¥INUS  1  $D 


~r~T — i     i     i — iiii — i — i — i — i     i     i — r— i — i     i 

0  5  10  15  20  25  30  35  40  45  50  55  60  05  70  75  80  65  00  05100 


TIME 


CASE  8  (2X2)  1  STABLE  ROOT 


Figure  4.4    Y  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  4.5    XI  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  4.6    X2  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  4.7    Yl  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  4.8    Y2  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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a)  The  M/W  trajectory  and  the  C/S  trajectory  were  very  similar  and  they  display 
the  same  trajectory  shape. 

b)  The  M/W  model  trajectory  was  outside  the  established  boundaries  for  the 
initial  portion  of  the  battle.  It  is  noted,  however,  that  the  area  bounded  by 
one  standard  deviation  is  extremely  small  in  this  region.  The  M/W  trajectory 
entered  the  bounded  area  as  the  two  curves  turned.  Both  models  reached  an 
equilibrium  point  in  the  vicinity  of  (200,200)  where  the  battle  would  end  in  a 
draw. 

c)  Examination  of  the  individual  unit  force  level  trajectories  in  Figures  E.5  and 
E.6  shows  that  the  M/W  trajectory  was  within  the  established  bounds  for  each 
force  as  a  function  of  time. 

3.  Case  Three 

Case  three  was  designed  to  test  the  lxl  scenario  for  the  mode  of  behavior 
where  one  unstable  equilibrium  exists  in  the  vicinity  of  (200,200).  The  following 
observations  are  noted  from  the  results  shown  in  Figures  E.7  through  E.9: 

a)  The  force  level  trajectories  of  both  models  displayed  the  expected  behavior  for 
a  case  where  one  unstable  equilibrium  exists.  Each  of  the  models  trajectories 
approached  the  equilibrium  and  then  turn  away  from  it. 

b)  Once  again  the  M/W  trajectory  was  well  within  one  standard  deviation  of  the 
C/S  trajectory. 

4.  Case  Four 

Case  four  was  designed  to  test  a  lxl  scenario  where  the  mode  of  behavior 
involves  two  equilibria  in  the  vicinity  of  (200,200)(U)  and  (366,100)(S).  The  following 
observations  are  noted  from  the  results  of  this  case  shown  in  Figures  E.10  through 
E.12: 

a)  The  results  show  the  trajectories  of  both  models  display  the  expected  behavior 
approaching  the  unstable  equilibrium  and  then  turning  away. 

b)  The  M/W  model  was  always  well  within  the  one  standard  deviation 
boundaries  for  this  case. 

E.       (1X2)  CASE  COMPARISONS 

One  case  was  developed  to  test  an  unbalanced  force  level  scenario  involving  one 
type  X  force  and  two  type  Y  force.  The  purpose  of  this  case  was  to  ensure  that  both 
models  perform  properly  for  an  unbalanced  number  of  forces  engaging  in  the  battle.  A 
1x2  scenario  was  chosen  for  this  case  and  the  results  are  presented  next. 
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1.  Case  Five 

Case  five  was  the  one  unbalanced  scenario  developed  and  the  mode  of 
behavior  tested  involved  one  unstable  equilibrium  in  the  vicinity  of  (300,200).  The 
following  observations  were  noted  from  the  results  shown  in  Figures  E.13  through 
E.18: 

a)  Both  models  displayed  the  expected  behavior  for  an  unstable  equilibrium.  The 
trajectories  for  both  models  approached  the  equilibrium  point  and  then  were 
turned  away. 

b)  The  M/W  trajectory  was  well  within  the  one  standard  deviation  boundaries, 
and  had  the  same  shape  as  the  C/S  trajectory.  It  is  interesting  to  note  the 
relative  speed  at  which  the  two  trajectories  moved  away  from  the  unstable 
equilibrium  point.  The  C/S  model  moved  away  much  faster  than  the  M/W 
model. 

c)  The  plot  of  each  force  vs.  time,  shown  in  Figures  E.14  through  E.18,  show  the 
two  models  results  separating  over  time,  and  the  M/W  trajectory  leaving  the 
bounded  area  toward  the  end  of  the  battle. 

d)  This  case  involving  a  1x2  scenario  with  an  unstable  equilibrium  is  similar  to 
case  three  which  was  a  lxl  scenario  involving  one  unstable  equilibrium. 

F.       (2X2)  CASE  COMPARISONS 

The  next  step  in  the  experiment  was  to  develop  cases  which  involved  two  types  of 
forces  on  each  side  and  compare  the  models  results  for  the  four  modes  of  behavior 
discussed.  Cases  six  through  nine  are  2x2  case  scenarios.  Their  results  are  shown  in 
Figures  E.19  through  E.46.  The  results  for  each  of  these  cases  will  be  discussed 
individually. 

1.  Case  Six 

Case  six  is  a  2x2  force  level  scenario  designed  to  test  the  models  when  there 
are  no  equilibria  in  the  state  space.  The  results  for  case  six  are  shown  in  Figures  E.19 
through  E.25.   The  following  observations  are  noted  from  the  results: 

a)  The  resulting  trajectories  shown  in  the  Total  Force  Level  Trajectory  plot  in 
Figure  E.19  show  both  models  behaved  similarly.  In  this  particular  case  the  Y 
force  was  the  winning  force.  The  expected  behavior  for  no  equilibria  was 
observed  in  this  case. 
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b)  The  M/W  model  results  were  within  the  defined  one  standard  deviation 
boundaries  during  most  of  the  battle,  leaving  the  boundaries  only  toward  the 
end  of  the  engagement. 

c)  Further  comparison  of  the  two  models  can  be  done  using  the  results  plotted 
for  the  individual  forces  in  Figures  E.20  through  E.25.  From  these  results  it  is 
evident  that  the  X  force  results  are  the  ones  where  the  models  differed  the 
most.   The  trajectories  all  had  the  same  general  trends. 

d)  This  2x2  case  involving  no  equilibria  is  similar  to  case  one  which  was  a  lxl 
case  scenario  involving  no  equilibria. 

2.  Case  Seven 

Case  seven  was  designed  to  test  the  models  results  with  two  type  X  forces  and 
two  type  Y  forces  when  two  equilibria  exist.  The  results  of  the  two  models  are  shown 
in  Figures  E.26  through  E.32.   The  following  observations  are  noted  for  case  seven: 

a)  The  resulting  trajectories  for  both  models  behaved  in  the  expected  manner. 
They  approached  the  unstable  equilibrium  and  then  turned  away  and  headed 
off.  In  this  particular  case  the  X  force  is  winning,  but  the  results  would  be 
different  if  the  force  levels  were  such  that  the  initial  point  was  above  the 
separatrix. 

b)  The  M/W  trajectory  had  the  same  shape  as  the  C/S  trajectory.  The  M/W 
trajectory  shown  in  Figure  E.26  was  outside  the  established  boundaries  during 
the  initial  portion  of  the  battle,  but  was  within  the  boundaries  after  they 
turned  away  from  the  equilibrium.  However,  examination  of  the  individual 
units  force  level  trajectories  shown  in  Figures  E.27  through  E.32  shows  the 
M/W  model  within  the  boundaries  for  the  X  force  and  just  outside  the 
boundaries  for  the  Y  force. 

c)  This  2x2  case  involving  two  equilibria  is  similar  to  case  four  which  was  a  lxl 
scenario  involving  two  equilibria. 

3.  Case  Eight 

Case  eight  was  designed  to  test  the  two  models  results  for  a  force  scenario 
involving  two  type  X  forces  and  two  type  Y  forces  with  one  stable  equilibrium  present 
in  the  state  space.  The  results  for  case  eight  are  shown  in  Figure  E.33  through  E.39. 
The  following  observations  are  noted  for  case  eight: 

a)       The  resulting  trajectories  for  the  M/W  and  C/S  models  both  exhibited  the 

expected   behavior  with   one   stable   equilibrium  present.     They   both  were 
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attracted  to  the  stable  equilibrium  point  and  remained  in  its  vicinity.   The  C'S 
model  shows  more  variability  due  to  its  stochastic  process. 

b)  The  M/W  trajectory  had  the  same  shape  and  remained  within  the  established 
boundaries  for  the  entire  battle. 

c)  This  2x2  case  involving  one  stable  equilibrium  is  similar  to  case  two  which  was 
a  lxl  scenario  involving  one  stable  equilibrium. 

4.  Case  Nine 

Case  nine  was  designed  to  test  the  two  models  results  for  a  2x2  force  level 
scenario  when  one  unstable  equilibrium  exists.  The  results  for  case  nine  are  shown  in 
Figures  E.40  through  E.46.   The  following  observations  are  noted  for  case  nine: 

a)  The  resulting  trajectories  for  each  of  the  models  behaved  in  the  expected 
manner  for  an  unstable  equilibrium. 

b)  The  M/W  trajectory  had  the  same  general  shape,  approaching  and  turning 
away  from  the  equilibrium,  but  it  was  not  within  the  designated  boundaries 
during  most  of  the  battle.  Reviewing  the  results  of  the  individual  units  shows 
that  the  two  models  did  not  behave  similarly  as  a  function  of  time. 

c)  This  2x2  case  involving  one  unstable  equilibrium  is  similar  to  case  three  which 
is  a  lxl  scenario  involving  one  unstable  equilibrium. 

G.       (3X3)  CASE  COMPARISONS 

The  last  force  level  scenario  developed  for  this  experiment  involves  three  type  X 
forces  and  three  type  Y  forces.  One  case,  labeled  as  a  3x3  case,  was  developed  for  this 
experiment.  There  is  no  discussion  of  the  two  models  results  for  each  of  the  four 
modes  of  behavior  discussed  earlier  for  the  3x3  scenario.  Rather,  one  case  was  run  to 
test  the  models  results  for  the  larger  3x3  scenario. 

1.  Case  Ten 

Case  ten  is  the  one  3x3  case  developed  for  comparing  the  M/W  model  results 
with  the  C/S  model  results.  The  following  observations  are  made  from  the  results 
shown  in  Figures  E.47  through  E.55: 

a)  The  M/W  trajectory  was  well  within  the  established  boundaries  throughout 
the  battle  and  the  two  trajectories  shapes  were  almost  identical. 

b)  The  M/W  and  C/S  trajectories  for  each  of  the  individual  units  are  very  similar. 
In  several  of  the  cases  the  two  trajectories  are  overlapping  for  most  of  the 
battle. 
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H.      SUMMARY  OF  RESULTS 

The  results  of  the  numerical  experiment,  with  ten  cases  involving  lxl,  1x2,  2x2, 
and  3x3  force  level  scenarios,  have  been  presented  and  briefly  discussed  in  this  chapter. 
Some  general  observations  from  these  results  will  serve  as  a  summary  of  the 
experiment. 

First,  the  results  for  the  ten  cases  support  the  idea  that  a  deterministic  attrition 
model  is  often  a  good  approximation  for  the  mean  result  of  a  stochastic  attrition 
model.  From  the  results  shown  in  the  previous  section,  only  one  of  the  ten  cases 
displayed  force  level  trajectories  where  the  M/W  model  was  significantly  different  from 
the  C/S  model.  Five  of  the  ten  cases  showed  the  M/W  trajectories  were  within  the 
established  one  standard  deviation  boundary  for  the  entire  battle  (cases  one,  three, 
four,  eight,  and  ten).  Of  the  remaining  cases  where  the  deterministic  M/W  model 
results  were  outside  the  boundary  for  a  portion  of  the  battle  shown  in  the  Total  Force 
Level  Trajectory  plot,  it  was  noted  that  the  individual  unit  force  level  trajectories  were 
within  the  established  one  standard  deviation  boundary  in  all  but  one  case  (case  nine). 

Secondly,  the  results  from  the  experiment  show  that  both  models  force  level 
trajectories  displayed  the  expected  mode  of  behavior  for  the  number  and  type  of 
equilibrium  points  which  were  present  in  all  cases.  This  shows  strong  support  for  the 
research  done  by  Professors  Moose  and  Wozencraft. 

This  experiments  results  are  certainly  not  proof  that  the  deterministic  model  is 
always  a  good  approximation  for  the  mean  result  of  a  stochastic  model.  However,  the 
results  do  indicate  they  often  may  be  a  good  approximation. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  results  of  this  research  showed  that  the  M/W  trajectory  was  a  good 
approximation  for  the  C/S  trajectory  in  nine  of  ten  cases.  This  is  an  indication  that 
deterministic  attrition  models  may  be  good  approximations  for  the  mean  of  stochastic 
attrition  models.  It  has  also  shown  that  the  modes  of  behavior  for  both  models 
trajectories  are  very  similar.  The  combination  of  these  two  results  provides  strong 
support  that  the  M/W  model  is  a  good  approximation  for  the  C/S  model.  These 
results  provide  important  support  for  the  argument  that  a  deterministic  attrition  model 
can  be  used  to  model  the  combat  attrition  process  which  is  generally  agreed  upon  as  a 
stochastic  process. 

This  simulation  provided  a  test  of  the  Moose/Wozencraft  model.  The  M/W 
model  was  tested  utilizing  four  different  force  composition  scenarios,  for  each  of  the 
four  types  of  trajectory  behavior  modes.  The  results  showed  the  M/W  models  results 
were  good  approximations  for  the  individual  unit  level  and  the  aggregate  force  levels  of 
the  C/S  model.  This  Generalized  Lanchester  Equation  type  model  is  simple  to  use  and 
understand.  It  allows  the  modeler,  or  military  planner,  to  model  an  aggregate  force  of 
as  many  different  type  forces  as  required  to  model  the  real-world  battlefield.  The  M/W 
model  has  real  potential  for  future  use  by  military  analysts. 

B.  COMMAND  AND  CONTROL  APPLICATIONS 

A  key  decision  in  the  Command  and  Control  decision  process  is  what  force 
structure  is  required  to  engage  a  proposed  enemy.  The  problem  is  to  determine  the 
optimum  force  size  and  composition  required  to  engage  an  enemy  and  win.  Optimizing 
one's  force  is  a  very  complex  problem,  involving  the  right  mixture  of  force  types  and 
the  right  number  of  each  of  the  forces.  An  example  would  be  when  a  commander  may 
need  a  tank  heavy  unit  for  a  specific  battle  in  order  to  ensure  that  his  force  will  win. 

A  model  such  as  the  M/W  model  can  serve  as  a  very  useful  decision  aid  for  the 
military  planners  and  combat  commanders  to  assist  with  optimization  of  forces  on  the 
battlefield.  Given  an  accepted  set  of  attrition-rate  coefficients,  the  model  can  provide 
the  capability  of  testing  for  the  optimum  mix  and  force  size  required  to  fight  and  win  a 
battle.     The   use    of  a   model   of  this   type   by   maneuver   units   would   require   the 
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commanders  and  their  staffs  to  ask  the  intelligence  community  for  specific  information 
on  the  enemies  current  force  structure  and  resupply  capability.  This  information  could 
be  used  with  the  model  to  assist  planners  and  commanders  to  decide  whether  to 
engage,  and  the  optimum  force  compostion  if  the  decision  is  made  to  do  so.  If  not, 
they  may  chose  to  with-draw  and  fight  on  better  terms.  Models  of  this  type  would  be 
beneficial  to  military  planners  responsible  for  developing  OPLANS  and  CONPLANS 
also.  It  would  allow  them  to  hypothesize  the  expected  enemy  force  structures  and  test 
proposed  force  structures  to  engage  the  enemy. 

C.       RECOMMENDATIONS 

After  completing  this  project  it  is  apparent  that  a  great  deal  of  future  research 
can  be  done  with  these  two  models  to  provide  a  more  complete  analysis.  The  results 
from  this  initial  survey  have  identified  areas  for  further  testing  which  would  provide  a 
more  complete  understanding  of  the  dynamics  of  the  two  models.  The  areas  listed 
below  would  provide  good  topics  for  future  work  with  these  models: 

1)  Utilize  the  models  to  locate  the  equilibrium  points  and  the  separatrix  for  cases 
where  an  unstable  equilibrium  point  exits  in  the  state  space.  Once  they  have 
been  located,  observe  the  two  models  results  in  the  vicinity  of  the  equilibrium 
point  for  several  different  inital  force  levels.  The  separatrix  is  a  very 
interesting  factor  which  should  be  observed  closely. 

2)  A  detailed  analysis  of  case  nine  of  this  experiment  is  recommended.  Case  nine 
was  the  one  case  which  was  run  that  did  not  behave  similarly  to  the  rest  of  the 
cases.  At  this  time  it  is  not  known  why,  and  therefore  deserves  further 
research. 

3)  One  of  the  observations  made  from  this  experiment  involved  a  difference  in 
the  attrition-rate  of  a  force  over  time  for  the  two  models.  This  could  be  a 
result  of  the  normalization  factor  used  for  the  aimed  fire  attrition  in  the  M/W 
model.  This  factor  was  used  to  ensure  that  an  individual  force  was  not 
allowed  to  fall  below  a  zero  level.  This  factor  may  cause  the  two  models  to 
have  slightly  different  attrition-rates  over  time.  This  can  be  tested  by 
removing  the  factor  and  changing  the  models  computer  program  to  stop  if  one 
force  falls  to  a  zero  level. 

4)  Test  the  model(s)  against  approved  combat  simulations. 
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Further  research  in  this  area  may  lead  to  an  attrition  model  which  is  accepted  as 
a  useful  model  for  modern  combat.  This  could  provide  a  useful  Command  and  Control 
decision  aid  for  use  by  commanders  and  their  staffs  to  optimize  force  structures  on  the 
battlefield. 
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APPENDIX  A 
EXEC  FILE  FOR  C/S  PROGRAM 


&TRACE  ON 

GLOBAL  TXTLIB  VFORTLIB  NONIMSL  IMSLSP  CMSLIB 

GLOBAL  LOADLIB  VFLODLIB 

FILEDEF  4  TERM 

FILEDEF  5  DISK  INTAW  DATA  A 

FILEDEF  6  DISK  TAW  FT06F001  A 

FILEDEF  7  DISK  TAW   FT07F001  A 

FILEDEF  8  DISK  TXTY  DATA  A 

FILEDEF  9  DISK  X1Y1  DATA  A 

FILEDEF  10  DISK  X2Y2  DATA  A 

FILEDEF  11  DISK  X3Y3  DATA  A 

FILEDEF  12  DISK  STEP  DATA  A 

FILEDEF  13  DISK  STATTF  DATA  A 

FILEDEF  14  DISK  STATF1  DATA  A 

FILEDEF  15  DISK  STATF2  DATA  A 

FILEDEF  16  DISK  STATF3  DATA  A 

LOAD  TAW 

START  MAIN 

CP  MSG  0731P  JOB  TAW  COMPLETED 
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APPENDIX  B 
C/S  COMPUTER  PROGRAM 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


PROGRAM  LANCHESTER 
6/18/86 
WRITTEN  TO  TEST  THE  MOOSE/WOZENCRAFT  THEORY 


NTYPX,Y 

NXS,YS(I) 

NX,Y(I) 

aaa(i,j; 

BBB(I,J 
CCC(I,J 
DDD(I,J] 

NRRR(I) 
NSSS(I) 

X,YPOS(K) 

X,YSTAT(K) 

NQX,Y(K) 

PKHX,Y(I,J) 
X,YAF(I,J) 

RN(K) 
RX,Y(K) 
RXS,YS(K) 
RXT,YT(K) 

PKX,Y(I,J) 
KAX,Y(I,J) 
KHX,Y(I,J) 
PKSX,Y(I) 

NDHX,Y(I) 
NDAX,Y(I) 
NDSX,Y(I) 

NSTOP 
X,YLOC 


NUMBER  OF  TYPES  OF  COMBATANTS  <INPUT> 
INITIAL  NUMBER  OF  EACH  TYPE  OF  COMBATANT 
NUMBER  OF  EACH  TYPE  OF  COMBATANT  ALIVE 


ATTRITION:  X(I' 

ATTRITION:  Y(J' 

ATTRITION:  Y(I 

ATTRITION:  X(I 


LOSSES  TO  AREA  FIRE 
SHOOTS  X 

LOSSES  TO  AREA  FIRE 
SHOOTS  Y 


NUMBER  COMBATANT  X(I) 
NUMBER  COMBATANT  Y(I) 


RESUPPLIED/TIMESTEP 
RESUPPLIED/TIMESTEP 


STORED  POSITION  OF  EACH  COMBATANT 

STORED  STATUS  (LIVE  OR  DEAD)  OF  EACH  COMBATANT 

STORED  TYPE  OF  EACH  COMBATANT 

PK:  X  SHOOTS  Y  TYPE  (I, J)  INTERACTION  <INPUT> 
AREA  FIRE  LOSS  COEFFICIENTS  FOR  X  <INPUT> 

LISTS  OF  RANDOM  VARIANTS,  CONTROLS  WHO  SHOOTS 

LISTS  OF  RANDOM  VARIANTS,  CONTROLS  HITS 

LISTS  OF  RANDOM  VARIANTS,  CONTROLS  SELF-LOSSES 

LISTS  OF  RANDOM  VARIANTS,  CONTROLS  TARGET  CHOICE 

PK:  X(I)  KILLS  Y(J)  THIS  TIMESTEP 
NR  OF  Y(J)  HIT  BY  X(I)  THIS  TIMESTEP,  AREA  FIRE 
NR  OF  Y(J)  HIT  BY  X(I)  THIS  TIMESTEP,  AIMED  FIRE 
PROBABILITY  OF  SELF-LOSS 

X,Y(I)  LOSSES  TO  AIMED  FIRE  THIS  TIMESTEP 
X,Y(I)  LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 
X,Y(I)  SELF  LOSSES  THIS  TIMESTEP 

MAXIMUM  NUMBER  OF  TIMESTEPS  <INPUT> 
POSSIBLE  POSITIONS  OF  COMBATANTS  <INPUT> 


IMPLICIT  REAL*8  (D) 

CHARACTER*40  IHEAD 

CHARACTERS   IANS ,  IYES 

DIMENSION  XPOS(3222) ,XSTAT(3222) ,NQX(3222) ,NX(3 

DIMENSION  YPOS(3222) ,YSTAT(3222) ,NQY(3222) ,NY(3 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


3),NRRR(3) 
),NSSS(3) 


XLOC(3),PKSX 

YLOC^KPKSYn, 

XAF(3,3),YAF(3,3 

PKX(3,3) .PKY(3,3 

PKHX ( 3 , 3 ) , PKHY ( 3 , 3 ) , KHX ( 3 , 3 ) , KHY ( 3 , 3 ) 

PKAX(3,3) ,PKAY(3,3),KAX(3,3),KAY(3,3) 

NDX(3) .NDY(3) ,NXS(3) ,NYS(3) 

NDHX(3),NDHY(3) 

NDAX(3),NDAY(3) 

NDSX(3) .NDSY(3) 

AAA(3,3),AAAAV(3,3),SAAAAV(3,3) 
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DIMENSION 
DIMENSION 
DIMENSION 

DIMENSION 
DIMENSION 
DIMENSION 

DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


,BBBAV 
,CCCAV 
, DDDAV 


SBBBAV 
SCCCAV 
SDDDAV 


RX(4222) .RY(4222) ,RN(9888) 

RXS(4222),RYS(4222) 

RXT(4222),RYT(4222) 


SNTX(212 
SNTY(212 


,SNX1(212),SNX2(212),SNX3(212 
SNY1(212) ,SNY2(212) ,SNY3(212 


SSNTX(212) ,SSNX1(212) ,SSNX2(212) ,SSNX3(212) 
SSNTY(212) ,SSNY1(212) ,SSNY2(212) ,SSNY3(212) 


ANTX(212) 
ANTY(212) 


SDNTX 
SDNTY 
PSNTX 
PSNTY 
TSNTX 
TSNTY 
CSTEP 


212 
'212 
'212' 

212 
'212 
'212' 

100' 


,ANX1(212) ,ANX2(212) ,ANX3(212) 
ANY1(212) ,ANY2(212) ,ANY3(212) 


,SDNX1 
,SDNY1 
,PSNX1 
,PSNY1 
,TSNX1 
,TSNY1 


212 
'212' 
'212 
'212' 
'212 
'212' 


,SDNX2i 
,SDNY2> 
,PSNX2' 
,PSNY2' 
,TSNX2' 
,TSNY2< 


212' 

212' 

212' 

212 

212' 

212' 


,SDNX3 
, SDNY3 
,PSNX3 
,PSNY3 
,TSNX3 
,TSNY3 


212' 
212' 
'212' 
'212' 
'212' 
'212' 


C 
C 


DATA  IYES/'Y1/ 


C 
C 

C 


CALL 

READ 

WRITE 
9410  FORMAT 

WRITE 

READ 
9405  FORMAT 

IF(IANS.EQ 
9400  FORMAT (A40) 


) 


EXCMS( 'CLRSCRN1 
'5,9400)  IHEAD 

4,9410) 

1  DO  YOU  WANT  TO  CHANGE  THIS  HEADING?  (Y/N)  (A40) ■ ) 

4,9400)  IHEAD 

4,9405)  IANS 

Al) 

I YES)   READ  (4,9400)  IHEAD 


ASK  USER  FOR  THE  NUMBER  OF  REPLICATIONS  DESIRED 

WRITE(4,9490) 
9490  FORMATf'   HOW  MANY  REPLICAS  DO  DESIRE?  (MAX  100)') 

READ(4,9495)  L 
9495  FORMAT (14) 

4,9480) 
9480  FORMAT ('    DO  YOU  WANT  FILES  6  AND  7  PRINTED?  ') 

4,9401)  IANS 
9401  FORMAT(Al) 

6,9400" 

7,9400' 

8,9400' 

'9,9400 

10,9400 

11,9400 

12,9400' 

13,9400 

14,9400 

15,9400 

16,9400, 

5,9500) 

5,9500) 


WRITE 
FORMAT 

READ 
FORMAT 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

READ 

READ 

IF(IANS.E 

WRITE(6,9 

WRITE (6,9501)  NTYPX , NX , XLOC , PKSX , NRRR 

ENDIF 

READ  (5,9505)  ( (PKHX(I , J) , J=l , 3) , 1=1 ,3) 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9521) 

WRITE (6, 9520)  ( (PKHX(I , J) , J=l , 3) , 1=1 , 3) 

ENDIF 

READ  (5,9505)  ( (XAF(I , J) , J=l ,3) , 1=1 ,3) 

IF ( IANS . EQ . IYES ) THEN 


IHEAD 
IHEAD 
IHEAD 
IHEAD 

IHEAD 

IHEAD 

IHEAD 

IHEAD 

IHEAD 

IHEAD 

IHEAD 
NSTOP 

NTYPX , NX , XLOC , PKSX , NRRR 
). IYES) THEN 

;o2; 
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WRITE(6,9523) 

WRITE (6, 9520)  ( (XAF(I , J) , J=l , 3) , 1=1 , 3) 

ENDIF 

READ  (5,9500)  NTYPY,NY,YLOC,PKSY,NSSS 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9503) 

WRITE (6, 9501)  NTYPY,NY, YLOC,PKSY,NSSS 

ENDIF 

READ  (5,9505)  ( (PKHY(I , J) , J=l ,3) , 1=1 ,3) 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9522) 

WRITE (6, 9520)  ( (PKHY(I , J) , J=l , 3) , 1=1 , 3) 

ENDIF 

READ  (5,9505)  ( (YAF(I , J) , J=l ,3) , 1=1 , 3) 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9524) 

WRITE (6, 9520)  ( (YAF(I , J) , J=l , 3) , 1=1 , 3) 

ENDIF 

9500  F0RMAT(4I5,3F5.1,3F7.3,3I5) 

9501  F0RMAT(4I5.3F5.1,3F7.3,3I5) 

9502  FORMAT (  '  #  X  TYPES  NUMBER  EACH    LOCATIONS      SELF-LOSS 
X        '    RE-SUPPLY     '  ) 

9503  FORMAT (  '  #  Y  TYPES  NUMBER  EACH    LOCATIONS      SELF-LOSS 
X        ■    RE-SUPPLY     ' ) 

9505  FORMAT (3F8. 5) 
9520  FORMAT (9F8. 5,/) 


9521  F0RMAT(/,'  KILL  PROBABILITIES  FOR  FORCE  X(I)  SHOOTS  Y(J)1,//, 

PK(1,1)  PK(1,2)  PK(1,3)  PK(2,1)  PK(2,2)  PK(2,3)\ 

PK(3,1)  PK(3,2)  PK(3,3)') 
KILL  PROBABILITIES  FOR  FORCE  Y(I)  SHOOTS  X(J)',//, 

PK(1,1)  PK(1,2)  PK(1,3)  PK(2,1)  PK(2,2)  PK(2,3)\ 

PK(3,1)  PK(3,2)  PK(3,3) ' ) 
AREA  FIRE  COEFFICIENTS  FORCE  Y(I)  SHOOTS  X(J)  ',//, 

AA(1,1)  AA(1,2)  AA(1,3)  AA(2,1)  AA(2,2)  AA(2,3)\ 

AA(3,1)  AA(3,2)  AA(3,3)') 
AREA  FIRE  COEFFICIENTS  FORCE  X(I)  SHOOTS  Y(J)  ',//, 

CC(1,1)  CC(1,2)  CC(1,3)  CC(2,1)  CC(2,2)  CC(2,3)'# 

CC  3,1)  CC  3,2)  CC  3,3  ') 


X 
X 

9522  FORMAT(/, ' 
X 
X 

9524  FORMAT(/, ' 
X 
X 

9523  FORMAT (/, ' 
X 
X 

c 

C  INITIALIZE  THE  RANDOM  NUMBER  GENERATORS 

C 

READ(5,9510)  DX,DY,DN 

DSEEDX=DX 

DSEEDY=DY 

DSEEDN=DN 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9511) 

WRITE (6, 9510)  DX,DY,DN,DSEEDX,DSEEDY,DSEEDN 

ENDIF 

9510  F0RMAT(3D15.8,3I11) 

9511  F0RMAT('     DOUBLE  PRECISION  RANDOM  VARIANT  SEEDS 
X       ■     INTEGER  SEED  VALUES ' ) 

C 

C     INITILAIZE  THE  STATISTICS  VARIABLES 

C 

DO  10  1=1,  212 

SNTX(I)=0 

SNX1(I)=0 

SNX2(I)=0 

SNX3(I)=0 

SNTY(I)=0 

SNYl(l)=0 

SNY2(I)=0 

SNY3(I)=0 

SSNTX(I)=0 

SSNX1(I)=0 

SSNX2(I)=0 

SSNX3(I)=0 

SSNTY(I)=0 
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SSNY1(I)=0 

SSNY2(I)=0 

SSNY3(I)=0 

ANTX(I)=0 

ANX1(I)=0 

ANX2(l)=0 

ANX3(I)=0 

ANTY(I)=0 

ANY1(I)=0 

ANY2(I)=0 

ANY3(I)=0 

SDNTX(I)=0 

SDNX1(I)=0 

SDNX2(I)=0 

SDNX3(l)=0 

SDNTY(I)=0 

SDNY1(I)=0 

SDNY2(I)=0 

SDNY3(I)=0 

PSNTX(I)=0 

PSNX1(I)=0 

PSNX2(I)=0 

PSNX3(I)=0  ■ 

PSNTY(I)=0 

PSNY1(I)=0 

PSNY2(I)=0 

PSNY3(I)=0 

TSNTX(I)=0 

TSNX1(I)=0 

TSNX2(I)=0 

TSNX3(I)=0 

TSNTY(I)=0 

TSNY1(I)=0 

TSNY2(I)=0 

TSNY3(I)=0 
10    CONTINUE 

DO  15  1=1,3 
NYS(I)=NY(I) 
NXS(I)=NX(I) 
15  CONTINUE 
MINST=99 
AVSTEP=0.0 
SDSTEP=0.0 
SSTEP=0.0 
SSSTEP=0.0 
DO  17  1=1,100 

CSTEP(I)=0.0 
17    CONTINUE 
C 

C     BEGIN  REPLICATIONS 
C 

DO  4000  N=1,L 

NSTEP=0 

NTX=0 

NTY=0 

DO  20  1=1,3 

NY(I)=NYS(I) 

NX(I)=NXS(I) 

IF(I.LE.NTYPX)   THEN 

NTX=NTX+NX(I) 
END  IF 
IF(I.LE.NTYPY)   THEN 

NTY=NTY+NY(I) 
ENDIF 
20  CONTINUE 
NTXS=NTX 
NTYS=NTY 
NXSTOP=50 
NYSTOP=50 
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c 

C  INITIALIZE  THE  TWO  FORCES 

C 

IMAX=0 

IMIN=1 

IP=0 

DO  200  J=1,NTYPX 

IMAX=IMAX+NX(J) 

DO  100  I=IMIN,IMAX 

IP=IP+1 

XPOS(IP)=XLOC(J) 

XSTAT(IP)=1.0 

NQX(IP)=J 
100  CONTINUE 

IMIN=IMAX+1 
200  CONTINUE 
C 

IMAX=0 

IMIN=1 

IP=0 

DO  400  J=1,NTYPY 

IMAX=IMAX+NY(J) 

DO  300  I=IMIN,IMAX 

IP=IP+1 

YPOS(IP)=YLOC(J) 

YSTAT(IP)=1.0 

NQY(IP)=J 
300  CONTINUE 

IMIN=IMAX+1 
400  CONTINUE 

DO  500  1=1,3 

NDX(I)=0 

NDY(I)=0 
500  CONTINUE 

DO  700  1=1,3 

DO  700  J=l,3 

AAAAV(I,J)=0.0 

BBBAV(I,J)=0.0 

CCCAV(I,J)=0.0 

DDDAV(I,J)=0.0 
700  CONTINUE 
C 

C  SETUP  A  FILE  FOR  FUTURE  PLOTTING 

C 

IF(IANS.EQ.IYES)THEN 

WRITE(7,9700) 
9700  FORMAT ('   #  NTXS  NTX  NX1   NX2  NX3  DX1  DX2  DX3 ' , 

X  ■  NTYS  NTY  NY1   NY2  NY3  DY1  DY2  DY3 ' ) 

WRITE (7, 97 10)  NSTEP ,NTXS ,NTX,NX,NDX,NTYS ,NTY,NY,NDY 

9710  FORMAT('+' ,14,2(515,314),/) 

WRITE(6,9711) 

9711  FORMAT ('  #   NTXS  NTX  NX1   NX2  NX3  RX1  RX2  RX3 ' , 
X  '  NTYS  NTY  NY1   NY2  NY3  RY1  RY2  RY3 ' ) 

WRITE (6,9710)  NSTEP , NTXS , NTX , NX , NRRR , NTYS , NTY , NY , NSSS 
END  IF 

TOTAL  NUMBER  OF  LIVE  COMBATANTS 
NUMBER  OF  COMBATANT  WHO  IS  FIRING 
TYPE  OF  COMBATANT  FIRING 
COUNT  OF  TARGETS  ATTACKED 
NUMBER  OF  TARGET  UNDER  FIRE 
TYPE  OF  TARGET  UNDER  FIRE 

THIS  IS  THE  BEGINNING  OF  THE  COMBAT  CYCLE 

1000  CONTINUE 

NSTEP=NSTEP+1 

NTT=NTX+NTY 

NGT=1.5*REAL(NTT) 
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C 

C 

NTT 

c 

NSX,Y 

c 

IX, Y 

c 

NTARX,Y 

c 

NCX,Y 

c 

JX,Y 

c 

c 

TH 

c 

NGX=1.3*REAL(NTX) 

NGY=1.3*REAL(NTY) 
C 
C  MAKE  THE  LISTS  OF  RANDOM  NUMBERS 


C 


CALL  GGUBS(DSEEDN,NGT,RN) 


CALL  GGUBS(DSEEDX,NGX,RX; 

CALL  GGUBS(DSEEDX,NGX,RXS) 

CALL  GGUBS(DSEEDX,NGX,RXT) 

CALL  GGUBS(DSEEDY,NGY,RY) 

CALL  GGUBS(DSEEDY,NGY,RYS) 

CALL  GGUBS(DSEEDY,NGY,RYT) 
C 

C  ZERO  STORAGE  FOR  THIS  TIMESTEP 

C 

NSX=0 

NTARX=0 

NSY=0 

NTARY=0 

KSTOPX=0 

KSTOPY=0 

DO  1040  1=1,3 

NDX(I)=0 

NDY(I)=0 

NDHX(I)=0 

NDHY(I)=0 

NDAX(l)=0 

NDAY(I)=0 

NDSX(l)=0 

NDSY(I)=0 

DO  1020  J=l,3 

KAX(I,J)=0 

KAY(I,J)=0 

KHX(I,J)=0 

KHY(I,J)=0 
1020  CONTINUE 
1040  CONTINUE 

DO  1060  1=1,3 

DO  1060  J=l,3 

PKAX ( I ,  J ) =XAF ( I , J) *NTY 

PKAY(I,J)=YAF(I,J)*NTX 
1060  CONTINUE 
1080  CONTINUE 

DO  1090  J=l,3 

DO  1090  1=1,3 

PKX(I,J)=PKAX(I,J)+PKHX(I,J) 

PKY(I,J)=PKAY(I,J)+PKHY(I,J) 
1090  CONTINUE 
C 

C  DO  AIMED  AND  AREA  FIRE  INTERACTIONS  FOR  THIS  TIMESTEP 

C 

C  CHOOSE  THE  COMBATANT  WHO  SHOOTS,  USING  RN(I) 

C 

C  MAKE  COMBATANT  NUMBER  (COUNT),  NSX,Y 

C  COMBATANT  MUST  BE  ALIVE 

C  DETERMINE  COMBATANT  TYPE,  LX,Y 

C  COMBATANT  MUST  NOT  BE  A  SELF-LOSS 

C 

C  INCREASE  THE  TARGET  COUNT,  NTARX,Y 

C  MAKE  NUMBER  OF  TARGET,  NCX,Y,  WITH  RXT,YT 

C  DETERMINE  TARGET  TYPE,  LX,Y 

C 

C  COMPARE:  PKX,Y(LX,LY)  >  RX,Y,  TARGET  IS  HIT 

C  PKX,Y(LX,LY)  <  RX,Y,  MISS 

C  PKHX,Y(LX,LY)  >  RX,Y,  HIT  BY  AIMED  FIRE 

C  PKHX,Y(LX,LY)  <  RX,Y,  HIT  BY  AIMED  FIRE 

C 
1100  CONTINUE 

DO  1180  M=1,NGT 
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c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 


IF(RN(M).LE.(REAL(NTX)/NTT))  THEN 
X  SHOOTS 


NSX=NSX+1 

IF(NSX.GT.NTX) 

IF(XSTAT(NSX).EQ.O.O) 

LX=NQX(NSX) 

IF(PKSX(LX).LT.RXS(NSX))  THEN 

NTARY=NTARY+1 

NCY=RYT(NTARY)*NTY+1 

LY=NQY(NCY) 

IF(PKX(LX,LY).GE.RX(NSX)) 

KILL  TARGET  NUMBER  NCY 


GO  TO  1120 
GO  TO  1180 


THEN 


IF(YSTAT(NCY) .GT.0.0)  THEN 
YSTAT(NCY)=0.0 

IF(PKHX(LX,LY).GE.RX(NSX))   THEN 
NDHY ( LY ) =NDHY ( L Y ) + 1 
KHX(LX,LY)=KHX(LX,LY)+1 
ELSE 

NDAY(LY)=NDAY(LY)+1 
KAX(LX,LY)=KAX(LX,LY)+1 
END  IF 
END  IF 
ENDIF 

GO  TO  1180 
ELSE 

XSTAT(NSX)=0.0 
NDSX(LX)=NDSX(LX)+1 
GO  TO  1180 
ENDIF 


ELSE 


Y  SHOOTS 


NSY=NSY+1 

IF(NSY.GT.NTY) 

IF(YSTAT(NSY).EQ.O.O) 

LY=NQY(NSY) 

IF(PKSY(LY).LT.RYS(NSY)) 

NTARX=NTARX+1 

NCX=RXT (NTARX) *NTX+1 

LX=NQX(NCX) 

IF(PKY(LY,LX).GE.RY(NSY)) 

KILL  TARGET  NUMBER  NCX 


GO  TO  1150 
GO  TO  1180 

THEN 


THEN 


THEN 


IF(XSTAT(NCX).GT.O.O)  THEN 
XSTAT(NCX)=0.0 
IF(PKHY(LY,LX).GE.RY(NSY)) 
NDHX ( LX ) =NDHX ( LX ) + 1 
KHY(LY,LX)=KHY(LY,LX)+1 
ELSE 

NDAX(LX)=NDAX(LX)+1 
KAY(LY,LX)=KAY(LY,LX)+1 
ENDIF 
ENDIF 
ENDIF 

GO  TO  1180 
ELSE 

YSTAT(NSY)=0.0 
NDSY(LY)=NDSY(LY)+1 
GO  TO  1180 
ENDIF 
ENDIF 


ALL  TYPE  X  COMBATANTS  HAVE  FINISHED  FIRING 
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c 

1120  CONTINUE 
KSTOPX=l 
GO  TO  1170 
C 

C  ALL  TYPE  Y  COMBATANTS  HAVE  FINISHED  FIRING 

C 
1150  CONTINUE 
KSTOPY=l 
C 

C  EVERY  ONE  HAS  FINISHED  FIRING 

C 
1170  CONTINUE 

IF(KSTOPX.EQ.l  .AND.  KSTOPY.EQ.l)      GO  TO  1200 
1180  CONTINUE 
C 

C  CALCULATE  THE  KILL  RATIOS 

C 
1200  CONTINUE 

DO  1290  J=1,NTYPX 
DO  1250  I=1,NTYPY 

NUMD=NX(J)*NY(I) 

IF  (NUND.GE.l)  THEN 

DDD(J,I)=REAL(KHX(J,I)*NTY)/REAL(NUMD) 
ELSE 

DDD(J,I)=0.0 
END  IF 

DDDAV(J,I)=DDDAV(J,I)+DDD(J,I) 
NUMB=NY(I)*NX(J) 
IF  (NUMB.GE.l)  THEN 

BBB(I,J)=REAL(KHY(I,J)*NTX)/REAL(NUMB) 
ELSE 

BBB(I,J)=0.0 
ENDIF 

BBBAV(I,J)=BBBAV(I,J)+BBB(I,J) 
NUMC=NX(J)*NY(I) 
IF  (NUMC.GE.l)  THEN 

CCC(J,I)=REAL(KAX(J,I))/REAL(NUMC) 
ELSE 

CCC(J,I)=0.0 
ENDIF 

CCCAV(J,I)=CCCAV(J,I)+CCC(J,I) 
NUMA=NY(I)*NX(J) 
IF  (NUMA.GE.l)  THEN 

AAA(I,J)=REAL(KAY(I,J))/REAL(NUMA) 
ELSE 

AAA(I,J)=0.0 
ENDIF 

AAAAV(I,J)=AAAAV(I,J)+AAA(I,J) 
1250  CONTINUE 
C 

C  TOTAL  KILLED  BY  ALL  MECHANISMS 

C 

NDX(J)=NDHX(J)+NDAX(J)+NDSX(J) 
NDY(J)=NDHY(J)+NDAY(J)+NDSY(J) 
1290  CONTINUE 

IF(IANS.EQ.IYES)THEN 

WRITE ( 6 , 9660 )  NSTEP , ( (KHX( I , J ) , J=l , 3 ) , 1=1 , 3 ) 
ENDIF 
9660  FORMAT(/,'  TIMESTEP',  14,/, 

X    10X,r  AIMED  FIRE  KILLS  OF  Y"S  BY  X  FORCE' 
X    ,//,'   D(l,l)   D(l,2)   D(l,3)   D(2,l)   D(2,2)   D(2,3)', 
X       '   D(3,l)   D(3,2)   D(3,3)',/,1X, 
X  918  /^ 

IF(IANS.EQ.IYESl)THEN 

WRITE (6, 9665)  ( (KHY(I , J) , J=l ,3) ,1=1 ,3) 
ENDIF 
9665  FORMAT(10X,'  AIMED  FIRE  KILLS  OF  X"S  BY  Y  FORCE  '  ,- 

X        //,'   B(l,l)   B(l,2)   B(l,3)   B(2,l)   B(2,2)   B(2,3)', 
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c 
c 
c 


X  '   B(3,l)   B(3,2)   B(3,3)\/,1X, 

X  918  / ) 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9662)  ( (KAX(I , J) , J=l ,3) , 1=1 ,3) 

ENDIF 
9662  FORMAT(10X  '  AREA  FIRE  KILLS  OF  Y"S   BY  X  FORCE  ' 

X        //,'   0(1,1   0(1,2)   C(l,3)   C(2,l)   C(2,2)   0(2,3)', 
X  '   0(3,1)   0(3,2)   0(3,3)',/, IX, 

X   /  918,/) 

IF(IANS.EQ.IYES)THEN 

WRITE (6, 9668)  ( (KAY(I , J) , J=l ,3) , 1=1 ,3) 

ENDIF 
9668  FORMAT(10X, '  AREA  FIRE  KILLS  OF  X"S   BY  Y  FORCE  ', 

X        //,'   A(l,l   A(l,2)   A  1,3)   A(2,l)   A(2,2)   A(2,3)', 
X  '   A(3,l)   A(3,2)   A(3,3)',/,1X, 

X  918,/) 

IF(IANS.EQ.IYES)THEN 

WRITE(6,9700) 

WRITE (6, 9710)  NSTEP,NTXS ,NTX,NX,NDX,NTYS ,NTY,NY,NDY 

WRITE(7,9651) 

ENDIF 
9651  FORNATC   D,B,C,A  BY  ROWS  IN  DESCENDING  ORDER  ') 

IF(IANS.EQ.IYES)THEN 

WRITE (7, 9650)  NSTEP, ( (KHX(I , J) , J=l , 3) , 1=1 ,3' 

WRITE (7, 9650)  NSTEP, ( (KHY(I , J) , J=l ,3) , 1=1 ,3' 

WRITE (7, 9650)  NSTEP, ( (KAX(I , J) , J=l , 3) , 1=1 ,3 

WRITE (7,9650)  NSTEP , ( (KAY( I , J ) , J=l , 3 ) , 1=1 , 3 

ENDIF 
9650  FORMAT (14, 918) 


1500 


1510 


CALCULATE  TOTAL  LOSSES  FOR  ALL  INTERACTIONS 

CONTINUE 

NDTX=0 

NDTY=0 

DO  1510  J=l,3 

NDTX=NDTX+NDX(j; 

NDTY=NDTY+NDY(J; 
CONTINUE 

I F ( I ANS . EQ . I YES ) THEN 

WRITE (6, 9670)  NDTX,NRRR,NDHX,NDAX,NDSX, 
X  NDTY,NSSS,NDHY,NDAY,NDSY 

WRITE (7, 9670)  NDTX,NRRR,NDHX,NDAX,NDSX, 


9670 


C 
C 
C 


X 

FORMAT ( 
X 
X 
X 
X 
X 
X 
X 
ENDIF 


NDTY,NSSS,NDHY,NDAY,NDSY 
TOTAL  X  FORCE  LOSS  WAS ',14  '   RESUPPLY  WAS:', 314,/, 
AIMED  FIRE  LOSSES: ' ,314,/, 
AREA  FIRE  LOSSES 
SELF  LOSSES 
TOTAL  Y  FORCE  LOSS  WAS ',14,' 
AIMED  FIRE  LOSSES:  '  ,314,/', 
AREA  FIRE  LOSSES 
SELF  LOSSES 


',314,/, 
',314,/, 


RESUPPLY  WAS:' ,314,/, 


,314./, 
\3I4J 


ADD  IN  RESUPPLIED  UNITS 


1600  CONTINUE 

DO  1650  K=l,3 

IF(K.LE.NTYPX  .AND.  NRRR(K) .NE .0)   THEN 
NADDX=NRRR(K) 
NX(K)=NX(K)+NADDX 
DO  1620  I=1,NADDX 
NTX=NTX+1 
XPOS(NTX)=XLOC(K) 
NQX(NTX)=K 
XSTAT(NTX)=1.0 
1620    CONTINUE 
ENDIF 
IF(K.LE.NTYPY  .AND.  NSSS(K) .NE .0)   THEN 
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c 
c 
c 


NADDY=NSSS(K) 

NY(K)=NY(K)+NADDY 

DO  1640  I=1,NADDY 

NTY=NTY+1 

YPOS(NTY)=YLOC(K) 

NOY(NTY)=K 

YSTAT(NTY)=1.0 

1640    CONTINUE 
END  IF 

1650  CONTINUE 


C 
C 

C 


C 

c 
c 


CLOSE  OUT  DESTROYED  UNITS 

1700  CONTINUE 

DO  1710  J=l,3 
NX(J)=0 
NY(J)=0 
1710  CONTINUE 
IX=0 

DO  1740  I=1,NTX 
IF(XSTAT(I) .NE.0.0)  THEN 

IX=IX+1 

NQX(IX)=NQX(I) 

XPOS(IX)=XPOS(I) 

XSTAT(IX)=XSTAT(I) 

K=NQX(I) 

NX(R)=NX(K)+1 
ENDIF 
1740  CONTINUE 
NTX=IX 
IY=0 

DO  1780  1=1, NTY 
IF(YSTAT(I).NE.0.0)  THEN 

IY=IY+1 

NQY(IY)=NQY(I) 

YPOS(IY)=YPOS(I) 

YSTAT(IY)=YSTAT(I) 

K=NQY(I) 

NY(K)=NY(K)+1 
ENDIF 
1780  CONTINUE 
NTY=IY 


PRINT  STATUS  FOR  EACH  TIME  STEP  IF  REQUESTED 

IF(IANS.EQ.IYES)THEN 

WRITE (7, 97 10)  NSTEP ,NTXS ,NTX,NX,NDX,NTYS ,NTY,NY,NDY 
10 


WRITE(6,97 
END  IF 


NSTEP , NTXS , NTX , NX , NRRR , NTYS , NTY , NY , NSSS 


ACCUMILATE  STATISTICAL  VARIABLES  FOR  EACH  TIME  STEP 


I=NSTEP 


SNTX 

SNX1 

SNX2 

SNX3 

SNTY 

SNY1 

SNY2 

SNY3 

SSNTX 

SSNX1 

SSNX2 

SSNX3 

SSNTY 

SSNY1 

SSNY2 

SSNY3 


=SNTXI 

=SNX1< 

=SNX2» 

=SNX3 1 

=SNTY< 

■-SNY1 

--SNY2 1 

=SNY3 1 

l=SSNTX< 

l=SSNXl' 

l=SSNX2i 

)=SSNX3i 

l=SSNTYl 

I=SSNY1 

l=SSNY2i 

l=SSNY3> 


I 

I 
I 

I' 

T 

I 

I' 

'I 


+NTX 

+NX(1' 

+NX(2' 

+NX(3' 

+NTY 

+NY(1' 

+NY(2' 

+NY(3' 


l+NTX**2 

l+NX(l)**2 

l+NX(2 

l+NX(3' 

i+NTY*' 

l+NY(l 

i+NY(2' 

t+NY(3' 
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c 
c 

c 

c 
c 
c 


c 
c 
c 


TERMINATE  IF  FORCES  ARE  ATTRITTED  BELOW  SPECIFIED  LEVEL 

IF(NTX.LE.NXSTOP  .OR.  NTY.LE.NYSTOP)   GO  TO  9120 

TERMINATE  ON  NUMBER  OF  CYCLES 

IF(NSTEP.LT.NSTOP)  GO  TO  1000 

9120  CONTINUE 

STEPS=REAL(NSTEP) 

CALCULATE  THE  AVERAGE  ATTRITION  COEFFICIENTS  PER  REPLICATION 

DO  9100  J=l,3 

DO  9100  1=1,3 

AAAAV(I,J)=AAAAV(I,J)/STEPS 

BBBAV(I,J)=BBBAV(I,J)/STEPS 

CCCAV(I,J)=CCCAV(I,J)/STEPS 

DDDAV(I , J)=DDDAV(I , J)/STEPS 

S AAAAV ( I , J ) =S AAAAV ( I , J ' 

=SBBBAV(I,J' 

=SCCCAV(I,J' 

=SDDDAV(I,J 


C 
C 
C 


SBBBAV(I 
SCCCAV(I 
SDDDAV(I 
9100  CONTINUE 

CSTEP(N)=NSTEP 


J 
J 
J 


+AAAAV(I,J 
+BBBAV(I,J 
+CCCAV(I,J 
+DDDAV(I,J, 


C 
C 
C 


ACCUMILATE  THE  STATISTICAL  VARIABLES  FOR  THE  NUMBER  OF  TIME  STEPS 

SSTEP=SSTEP+CSTEP (N) 
SSSTEP=SSSTEP+(CSTEP(N)**2) 

WRITE  THE  NUMBER  OF  REPLICATIONS  AND  THE  NUMBER  OF  STEPS  TO  FILE 


C 
C 
C 


WRITE (12, 8510 )N,CSTEP(N) 
8510  FORMAT ( 13, F4.0) 

KEEP  TRACK  OF  MINIMUM  NUMBER  OF  STEPS  FOR  ALL  REPLICATIONS 


C 
C 
C 


IF(MINST.GT.CSTEP(N))THEN 

MINST=CSTEP(N) 
ENDIF 
4000  CONTINUE 

CALCULATE  AVERAGE /STANDARD  DEV.OF  FORCE  LEVELS  FOR  EACH  TIME  STEP 


B=REAL(L 
RD=1.0/B 
RDD=1.0/ 
DO  4400 


C 
C 
C 
C 


ANTX 

ANX1 

ANX2 

ANX3 

ANTY 

ANY1 

ANY2 

ANY3 

SDNT 

SDNX1 

SDNX2 

SDNX3 

SDNTY 

SDNY1 

SDNY2 

SDNY3 


) 

(B*(B-1)) 
I=1,MINST 


=SNTX( 
=SNX1< 
=SNX2i 
=SNX3i 
=SNTYl 
=SNY1' 
=SNY2i 
=SNY3 i 


>*RD 
i*RD 
i*RD 
l*RD 
i*RD 
i*RD 
i*RD 
i*RD 


=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABSi 
=SQRT(ABS< 


B*SSNTXI 

'J) 

i-SNTXl 

'X 

i**2^ 

i*RDD)) 

B*SSNX1< 

I 

i-SNXll 

I 

i**2J 

i*RDD)) 

B*SSNX2| 

I 

I-SNX2I 

I 

i**2 

i*RDD)) 

B*SSNX3< 

I 

I-SNX3I 

I 

1**2 

i*RDD)) 

B*SSNTY< 

I 

i-SNTYl 

I 

1**2 

i*RDD)) 

B*SSNY1 ( 

I 

i-SNYl! 

I 

1**2 

i*RDD)) 

B*SSNY2< 

I 

•-SNY2I 

I 

1**2! 

l*RDD)) 

B*SSNY3l 

i] 

I-SNY3I 

[i; 

1**2] 

i*RDD)) 

CALCULATE  THE  VALUES  FOR  PLUS  AND  MINUS  ONE  STANDARD  DEVIATION 
FROM  THE  MEAN  VALUES  FOR  PLOTTING 
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c 
c 

c 


PSNTXl 

PSNX1' 

PSNX2i 

PSNX3I 

PSNTYl 

PSNY1 

PSNY2' 

PSNY3i 

TSNTXi 

TSNX1 

TSNX2' 

TSNX3i 

TSNTYi 

TSNY1 

TSNY2« 

TSNY3' 


=ANTXI 
=ANX1' 
=ANX2 I 
=ANX3 i 
=ANTYl 
=ANY1 
=ANY2 i 
=ANY3 I 
=ANTXl 
=ANX1 i 
=ANX2i 
=ANX3 i 
=ANTYi 
=ANY1 i 
=ANY2 I 


+SDNTXI 

+SDNX1 

+SDNX2' 

+SDNX3' 

+SDNTYI 

+SDMY1 I 

+SDNY2' 

+SDNY3< 


-(SDNTX(I)) 


SDNX1 
-SDNX2> 
-SDNX3i 
-SDNTYi 
-SDNY1 
-SDNY2 


9200 
C 
C 
C 


i=ANY3(I)-SDNY3i 

CALCULATE  THE  AVERAGE  ATTRITION  COEFFICIENTS  OVER  ALL  REPLICATIONS 

DO  9200  J=l,3 
DO  9200  K=l,3 

AAAAV ( K , J ) =S AAAAV ( K , J ) *RD 

BBBAV(K,J)=SBBBAV(K,J  *RD 

CCCAV(K,J)=SCCCAV(K; J)*RD 

DDDAV ( K , J ) =SDDDAV ( K , J ) *RD 
CONTINUE 

PRINT  STATISTICS  TO  FILES  8,9,10,11  FOR  PLOTTING 


C 
C 
C 


WRITE(8,8100)  I,ANTX(I) ,ANTY(I) ,PSNTX(I) ,PSNTY(I) ,TSNTX(I) , 
XTSNTY(I) 

WRITE(9,8100)  I,ANX1(I) ,ANY1(I) ,PSNX1(I) ,PSNY1(I) ,TSNX1(I) , 
XTSNYl(I) 

WRITE (10, 8100)  I, ANX2 ( I ) , ANY2 ( I ) , PSNX2 ( I ) , PSNY2 ( I ) , TSNX2 ( I ) , 
XTSNY2(I) 

WRITE (11,8100)  I , ANX3 ( I ) , ANY3 ( I ) , PSNX3 ( I ) , PSNY3 ( I ) , TSNX3 ( I ) , 
XTSNY3(I) 
8100  F0RMAT(I3,6(F8.2,1X)) 

WRITE  AVG'S  AND  STANDARD  DEVIATIONS  TO  FILE  FOR  TABLES 


C 
C 
C 


WRITE(13, 8150)1, ANTX(I) ,SDNTX(I ) , ANTY(I ) , SDNTY(I 
WRITE (14,8150)1, ANX1 ( I ) , SDNX1 ( I ) , ANY1 ( I ) , SDNY1 ( I 
WRITE (15,8150)1, ANX2 ( I ) , SDNX2 ( I ) , ANY2 ( I ) , SDNY2 ( I 
WRITE (16,8150)1, ANX3 ( I ) , SDNX3 ( I ) , ANY3 ( I ) , SDNY3 ( I 

8150  F0RMAT(I3,4(F8.2,1X)) 

4400  CONTINUE 


C 

c 

c 


CALCULTE  THE  AVG  AND  STANDARD  DEVIATION  FOR  THE  #  OF  NSTEPS  PER  RUN 

AVSTEP=SSTEP'lcRD 
SDSTEP=SQRT(ABS((B*SSSTEP-SSTEP**2)*RDD)) 
WRITE (12, 8500 )MINST , AVSTEP , SDSTEP 
8500  F0RMAT(I3,2F7.2) 

EXPLAIN  CAUSE  OF  TERMINATION 

9000  CONTINUE 

WRITE(7,9656) 
9656  FORMATS   A,B.C,D  BY  ROWS  IN  DESCENDING  ORDER  ') 

WRITE (7, 9655)  ( (AAAAV(I , J) , J=l , 3) , 1=1 , 3) 

WRITE (7, 9655)  ( (BBBAV(I , J) , J=l ,3) , 1=1 ,3) 
9655  FORMAT ('  AV  ',9F8.5) 


WRITE (6, 9695)  ( (AAAAV(I , J) , J=l ,3) , 1=1 ,3) 
-)  ((BBBAV(I,J),J=1,3),I=1,3) 


WRITE(6,9697 
9690  FORMAT(/,'  TlMESTEP ' ,  14,/ 

X  10X,'  AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 
X       ■  CC(1,1)  CC(1,2)  CC(1,3)  CC(2,1)  CC(2,2)  CC 
1  CC(3,1)  CC(3,2)  CC(3,3)',/, 


X 


(  X  E+03)',//, 
(2,3)', 
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9695 


9697 


9698 


9699 


X 

WRITE 
WRITE 
WRITE 
WRITE 

FORMAT 

X     ' 

X 

X 

X 

FORMAT 

X     ' 

X 

X 

X 

FORMAT 

X     ■ 

X 

X 

X 

FORMAT 

X     ' 

X 

X 

X 


,J),J=1,3),I=1,3 
,J),J=1,3),I=1,3 
,J),J=1,3),I=1,3 
,J)/J=l,3),I=l,3l 

COEFFICIENTS  FOR  Y  FORCE 
AA(1,3)  AA(2,1)  AA(2,2)  AA(2,3) 
AA(3,3)»,/,1X, 


3P,9F8.5,/) 
7,9655)  ((CCCAV(I,, 
7,9655)  ((DDDAV(I,, 

6.9698)  ((CCCAV(I,, 

6.9699)  ((DDDAVCI,, 
10X, 
AVERAGE  ATTRITION 

1  AA(1,1)  AA(1,2) 
1  AA(3,1)  AA(3,2) 
9F8.5,/) 
(10X, 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 
1  BB(1,1)  BB(1,2)  BB(1,3)  BB(2,1)  BB(2,2)  BB(2,3) 
'  BB(3,1   BB(3,2)  BB(3 , 3) ' , / , IX, 
9F8.5,/) 
(10X, 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 
1  CC(1,1)  CC(1,2' 
1  CC(3,1)  CC(3,2' 
,     9F8.5,/) 
(10X, 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 
1  DD(1,1)  DD(1,2)  DD(1,3)  DD(2,1)  DD(2,2)  DD(2,3) 
1  DD(3,1   DD(3,2)  DD(3 , 3) ' , / , IX, 
9F8.5,/) 


'CC(1,3)  CC(2,1)  CC(2,2)  CC(2,3) 
CC(3,3)',/,1X, 


,//, 


,//, 


,//, 


,//, 


9620 
9630 
9640 


IF(NTX.LT.NXSTOP)   THEN 

WRITE  (6,9620)  NSTEP,NTX,NTXS 
WRITE  (7,9620)  NSTEP,NTX,NTXS 
ENDIF 
IF(NTY.LT.NYSTOP)   THEN 

WRITE  (6,9630)  NSTEP,NTY,NTYS 
WRITE  (7,9630)  NSTEP,NTY,NTYS 
ENDIF 

IF  (NSTEP  .EQ.  NSTOP)   THEN 
WRITE  (6,9640)  NSTOP 
WRITE  (7,9640)  NSTOP 
ENDIF 
FORMAT ('  STOPPED  IN  STEP  ',15,'  BY  FORCE  X1,/, 

X       '  CURRENT  X  STRENGTH  =',15,',  INITIAL  X  STRENGTH' , 15 ,//) 
FORMAT ('  STOPPED  IN  STEP  ',15,'  BY  FORCE  Y\/, 

X       ■  CURRENT  Y  STRENGTH  =',15,',  INITIAL  Y  STRENGTH' , 15 ,//) 
FORMATC  STOPPED  BY  PROGRAM  AT  NSTOP  =  ',15,//) 
END 
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APPENDIX  C 
EXEC  FILE  FOR  M/W  PROGRAM 


&TRACE  ON 

GLOBAL  TXTLIB  VFORTLIB  NONIMSL  IMSLSP  CMSLIB 

GLOBAL  LOADLIB  VFLODLIB 

FILEDEF  4  TERM 

FILEDEF  5  DISK  INMWM  DATA  A 

FILEDEF  6  DISK  MWM   STLIST  A 

FILEDEF  8  DISK  TOTAL  DATA  A 

FILEDEF  9  DISK  FORI  DATA  A 

FILEDEF  10  DISK  FOR2  DATA  A 

FILEDEF  11  DISK  FOR3  DATA  A 

LOAD  MWM 

START  MAIN 

CP  MSG  0731P  JOB  MWM  COMPLETED 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


APPENDIX  D 
M/W  COMPUTER  PROGRAM 

PROGRAM  WRITTEN  TO  RUN  THE  GENERALIZED  LANCHESTER 

EQUATION  MODEL  DESIGNED  BY 
PAUL  H.  MOOSE  AND  JACK  M.  WOZENCRAFT 

1/08/87 
WRITTEN  TO  TEST  THE  MOOSE/WOZENCRAFT  THEORY 


NTYPX 

NTYPY 

TNX(I) 

TNY(J) 

STNX,Y 

TTNX,Y 

AAA(J,I 
BBB(J,I' 

cccci^1 

DDD(I,j; 

SAAA(I 
SBBB(I' 
SCCC(J' 
SDDD(J' 


NUMBER  OF  TYPES  OF  X  COMBATANTS  <INPUT> 

NUMBER  OF  TYPES  OF  Y  COMBATANTS  <INPUT> 

NUMBER  OF  EACH  TYPE  OF  X  COMBATANTS 

NUMBER  OF  EACH  TYPE  OF  Y  COMBATANTS 

VARIBLE  TO  SUM  TOTAL  FORCE  LEVEL  PER  TIMESTEP 

TOTAL  FORCE  LEVEL  PER  TIMESTEP 

ATTRITION:  X(I)  AREA  FIRE  LOSSES  DUE  TO  Y(J) 
ATTRITION:  X(I)  AIMED  FIRE  LOSSES  DUE  TO  Y(J) 
ATTRITION:  Y(J)  ARERA  FIRE  LOSSES  DUE  TO  X(I) 
ATTRITION:  Y(J)  AIMED  FIRE  LOSSES  DUE  TO  X(I) 


ATTRITION:  X(I 

ATTRITION:  X(I 

ATTRITION:  Y(J' 

ATTRITION:  Y(J' 


SUM  OF  THE  AREA  LOSSES 
SUM  OF  THE  AIMED  FIRE  LOSSES 
SUM  OF  THE  AREA  FIRE  LOSSES 
SUM  OF  THE  AIMED  FIRE  LOSSES 


NRRR(I 
NSSS 


(5) 


NUMBER  COMBATANT  X(I) 
NUMBER  COMBATANT  Y(J) 


RESUPPLIED/TIMESTEP 
RESUPPLIED/TIMESTEP 


UUU(I) 
VW(J) 

AMLX(I 
AMLY(J 


SELF-ATTRITION  COEFFICIENT  FOR  X  FORCE 
SELF-ATTRITION  COEFFICIENT  FOR  Y  FORCE 


Ul 


LOSSES  TO  AIMED 
LOSSES  TO  AIMED 


FIRE  THIS 
FIRE  THIS 


TIMESTEP 
TIMESTEP 


ARLX 
ARLY 


BS 


SLX(I 
SLY  (J 


TNDX(I) 
TNDY(J) 

NSTOP 
XSTOP 
YSTOP 


m 
m 


LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 
LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 


SELF  LOSSES  THIS 
SELF  LOSSES  THIS 


TOTAL  CHANGE  IN  X(I' 
TOTAL  CHANGE  IN  Y(j] 


TIMESTEP 
TIMESTEP 

FORCE  LEVEL  THIS 
FORCE  LEVEL  THIS 


TIME 
TIME 


STEP 
STEP 


MAXIMUM  NUMBER  OF  TIMESTEPS  <INPUT> 
X  FORCE  LEVEL  DEFINED  AS  THE  BREAK  PT. 
Y  FORCE  LEVEL  DEFINED  AS  THE  BREAK  PT. 


IMPLICIT  REAL*8  (D) 

CHARACTER*40  IHEAD 
CHARACTER*1   IANS , IYES 

DIMENSION  TNX(3).TNY(3) 
DIMENSION  NRRR(3),NSSS(3' 
DIMENSION  AMLX(3),AMLY(3' 
DIMENSION  ARLX(3),ARLY(3' 
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c 

c 

c 
c 


9410 


9405 

9400 
9401 


9500 
9501 
9502 

9503 

9505 
9520 
9521 


9522 
9524 
9523 


C 
C 
C 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


SLX(3) ,SLY(3) 

TNDX(3),TNDY(3) 

AAA(3,3),SAAA(3 

BBB(3,3),SBBB 

CCC(3,3),SCCC 

DDD(3,3),SDDD 

UUU(3),WV(3) 


DATA  IYES/'Y1/ 


CALL 

READ 

WRITE 
FORMAT 

WRITE 

READ 
FORMAT 

IF(IANS.EQ 
FORMAT (A40) 
FORMAT 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

READ 

READ 

WRITE 

WRITE 

READ 

WRITE 

WRITE 

READ 

WRITE 

WRITE 

READ 

WRITE 

WRITE 

READ 

WRITE 

WRITE 

READ 

WRITE 

WRITE 
FORMAT 
FORMAT 
FORMAT 
X 

FORMAT ( 
X 


) 


EXCMS( 'CLRSCRN' 
5,9400)  IHEAD 
4,9410) 

1  DO  YOU  WANT  TO  CHANGE  THIS  HEADING?  (Y/N)  (A40) ' ) 
4,9400)  IHEAD 
4,9405)  IANS 
Al) 

IYES)   READ  (4,9400)  IHEAD 


IHEAD 
IHEAD 
IHEAD 

IHEAD 

IHEAD 
NSTOP 
NTYPX,TNX,UUU,NRRR 

NTYPX,TNX.UUU,NRRR 
((DDD(I,J),J=1,3),I=1,3) 


((DDD(I,J 
((CCC(I,J 


J=l 
J=l 


FORMAT 

FORMAT 

FORMAT (/, 

X 

X 

FORMAT(/, 

X 

X 

FORMAT(/, 

X 

X 

FORMAT (/, 

X 

X 


Al) 

6,9400) 
8,9400) 
9,9400) 
10,9400) 
11,9400) 
5,9500' 
5,9500 
6,9502 
6,9501 
5,9505 
6,9521 
6,9520 
5,9505 
6,9523 
6,9520 
5,9500 
6,9503 
6,9501 
5,9505 
6,9522 
6,9520 
5,9505 
6,9524 

6,9520)  ((AAA(I,J) 
I5,3F5.0,3F7.3,3I5 
I5,3F7.2,3F7.3,3I5 
1  #  X  TYPES  NUMBER 

RE-SUPPLY     ' ) 
#  Y  TYPES  NUMBER 

RE-SUPPLY     ' ) 


),I=1,3) 
),I=1,3) 


((CCC(I,J),J=1,3),I=1,3) 
NTYPY,TNY,VVV,NSSS 

NTYPY,TNY.VW,NSSS 
((BBB(I,J),J=1,3),I=1,3) 


((BBB(I, J' 


!aaa(i;j; 


),J=1,3),I=1,3) 
),J=1,3),I=1,3) 


J=l,3),l=l,3) 


EACH  SELF-LOSS 
EACH  SELF-LOSS 


'3F8.5) 
;9F8.5,/) 


KILL  PROBABILITIES  FOR  FORCE  X(l) 


DD(1,1)  DD(1,2)  DD(1,3)  DD(2 , 1)  DD(2 ,2)  DD( 
DD(3,1)  DD(3,2)  DD(3,3)') 


DD(3,2 

KILL  PROBABILITIES  FOR 
BB(1,1)  BB(1,2)  BB(1,3) 
BB(3,1)  BB(3,2)  BB(3,3)' ) 

AREA  FIRE  COEFFICIENTS  FORCE 


SHOOTS  Y(JV ,//, 


FORCE  Y(I) 

BB(2,1)  BB(2,2)  BB( 


SHOOTS  X(JV ,//, 


AA(1,1)  AA(1,2)  AA(1,3)  AA(2,1) 
AA(3,1)  AA(3,2)  AA(3,3)') 
AREA  FIRE  COEFFICIENTS  FORCE  X(I) 
CC(1,1)  CC(1,2)  CC(1,3)  CC(2,1) 
CC(3,1)  CC(3,2)  CC(3,3)') 


Y(I)  SHOOTS  X(J)  ' ,//, 
AA(2,2)  AA(2,3)' , 


SHOOTS  Y(J)  ' ,//, 
CC(2,2)  CC(2,3)' , 


INITIALIZE  THE  INTIAL  FORCES  LEVELS 
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INTX=0 

INTY=0 

DO  1100  1=1,3 

INTX=INTX+TNX(I) 
INTY=INTY+TNY(I) 
1100  CONTINUE 

TTNX=INTX 

TTNY=INTY 

XSTOP=50.0 

YSTOP=50.0 

MINST=99 
C 

C  SETUP  A  FILE  FOR  TRACING  EACH  STEPS  RESULTS 

C 

WRITE(6,9711) 

9711  FORMAT ('  #   INTX  TTNX  TNX1  TNX2  TNX3  RX1  RX2  RX3 ' , 
X  '  INTY  TTNY  TNY1  TNY2  TNY3  RY1  RY2  RY3 ■ ) 

WRITE (6, 97 12)  NSTEP , INTX, TTNX, TNX,NRRR, INTY, TTNY, TNY,NSSS 

9712  FORMAT ( 14, 2 ( 13, 4F7. 2, 313),/) 

WRITE(6,9700) 
9700  FORMAT ('   #  INTX  TTNX  TNX1  TNX2  TNX3 ■ , 

X  '  INTY  TTNY  TNY1  TNY2  TNY3 ' ) 

WRITE ( 6 , 97 1 0 )  NSTEP . INTX , TTNX , TNX , INTY , TTNY , TNY 

9710  FORMAT(I4,2(I5,4F7.2),/) 
C     END  IF 
C 
C 

C  THIS  IS  THE  BEGINNING  OF  THE  COMBAT  CYCLE 

C 

DO  4000  N=l,NSTOP 
NSTEP=N 
C 

C  ZERO  STORAGE  FOR  THIS  TIMESTEP 

C 

STNX=0.0 
STNY=0.0 

DO  1015  I=1,NTYPX 
SAAA(I)=0.0 
SBBB(I)=0.0 
1015   CONTINUE 

DO  1020  J=1,NTYPY 
SCCC(J)=0.0 
SDDD(J)=0.0 
1020   CONTINUE 
C 

C      CALCULATE  THE  LOSSES  DUE  TO  AREA  FIRE  FOR  THIS  TIMESTEP 
C 

DO  1110  I=1,NTYPX 
DO  1105  J=1,NTYPY 

SAAA(I)=SAAA(I)+(AAA(J,I)*TNY(J)) 
SCCC(J)=SCCC(J)+(CCC(I,J)*TNX(I)) 
1105   CONTINUE 
1110   CONTINUE 

DO  1115  I=1,NTYPX 

ARLX(I)=(TNX(I)*SAAA(I)) 
1115   CONTINUE 

DO  1120  J=1,NTYPY 

ARLY(J)=TNY(J)*SCCC(J) 
1120   CONTINUE 
C 

C    CALCULATE  THE  LOSSES  DUE  TO  AIMED  FIRE  PER  TIMESTEP 
C 

DO  1130  I=1,NTYPX 
DO  1125  J=1,NTYPY 

SBBB(I)=SBBB(I)+(BBB(J,I)*TNY(J)) 
SDDD(J)=SDDD(J)+(DDD(I,J)*TNX(l)) 
1125   CONTINUE 
1130   CONTINUE 

DO  1135  I=1,NTYPX 
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AMLX(I)=(TNX(I)/TTNX)*SBBB(I) 
1135   CONTINUE 

DO  1140  J=1,NTYPY 

AMLY(J)=(TNY(J)/TTNY)*SDDD(J) 
1140   CONTINUE 
C 

C    CALCULATE  THE  SELF  LOSSES  PER  TIMESTEP 
C 

DO  1142  I=1,NTYPX 

SLX(I)=TNX(I)*UUU(I) 
1142   CONTINUE 

DO  1145  J=1,NTYPY 

SLY(J)=TNY(J)*WV(J) 
1145   CONTINUE 
C 

C    CALCULATE  THE  TOTAL  KILLED  BY  ALL  MECHANISMS 
C 

DO  1150  I=1,NTYPX 

TNDX(I)=-SLX(I)-ARLX(I)-AMLX(I)+NRRR(I) 
1150   CONTINUE 

DO  1155  J=1,NTYPY 

TNDY(J)=-SLY(J)-ARLY(J)-AMLY(J)+NSSS(J) 
1155   CONTINUE 
C 

C    CALCULATE  THE  FORCE  LEVELS  AFTER  THIS  TIMESTEP 
C 

DO  1160  I=1,NTYPX 

TNX(I)=TNX(I)+TNDX(I) 
1160   CONTINUE 

DO  1165  J=1,NTYPY 

TNY(J)=TNY(J)+TNDY(J) 
1165   CONTINUE 
C 

C   CALCUTE  THE  TOTAL  X,Y  FORCE  LEVELS 
C 

DO  1170  I=1,NTYPX 
STNX=STNX+TNX(I) 
1170   CONTINUE 

TTNX=STNX 
DO  1175  J=1,NTYPY 
STNY=STNY+TNY(J) 
1175   CONTINUE 

TTNY=STNY 
C 
C 

C        PRINT  STATUS  FOR  EACH  TIME  STEP 
C 

WRITE (6, 67 10)  NSTEP. INTX,TTNX,TNX, INTY,TTNY,TNY 

6710  FORMAT(I4,2(I5,4F7.2),/) 
C 
C 
C 

C     PRINT  VALUES  TO  FILES  8,9,10,11  FOR  PLOTTING 
C 
C 
C 

WRITE (8, 8100)  N,TTNX,TTNY 

WRITE(9,8100)    N,TNX(1).TNY(1) 

WRITE(10,8100)   N,TNX(2),TNY(2; 

WRITE(11,8100)    N^NXUhTNYO' 
8100   F0RMAT(I3,2(F8.2,1X)) 
C 

C      TERMINATE  IF  FORCES  ARE  ATTRITTED  BELOW  SPECIFIED  LEVEL 
C 

IF(TTNX.LE.XSTOP  .OR.TTNY.LE .YSTOP)   GO  TO  5000 
C 

C   RETURN  AND  CALULATE  FORCE  LEVELS  FOR  ANOTHER  TIMESTEP 
C 
4000  CONTINUE 
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5000  CONTINUE 
C 

C  EXPLAIN  CAUSE  OF  TERMINATION 

C 
C 

IF(TTNX.LT.XSTOP)   THEN 

WRITE  (6,9620)  NSTEP,TTNX, INTX 
ENDIF 
IF(TTNY.LT.YSTOP)   THEN 

WRITE  (6,9630)  NSTEP,TTNY, INTY 
ENDIF 
IF  (MSTEP  .EQ.  NSTOP)   THEN 

WRITE  (6,9640)  NSTOP 
ENDIF 
9620  FORMAT ('  STOPPED  IN  STEP  ',15,'  BY  FORCE  X',/, 

X       '  CURRENT  X  STRENGTH  = ' , F5 . 2 , ' ,  INITIAL  X  STRENGTH' , 15 ,// ) 
9630  FORMAT ('  STOPPED  IN  STEP  ',15,'  BY  FORCE  Y',/, 

X       '  CURRENT  Y  STRENGTH  =',F5.2,',  INITIAL  Y  STRENGTH1  , 15  ,// ) 
9640  FORMAT (■  STOPPED  BY  PROGRAM  AT  NSTOP  =  ' ,15,//) 
END 
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APPENDIX  E 
FORCE  LEVEL  TRAJECTORIES 

This  Appendix  is  designed  to  provide  the  reader  with  a  complete  set  of  all  of  the 
results  of  this  experiment.  The  input  data  sets  for  each  model  are  presented  in  tables 
and  then  each  of  the  force  level  trajectories  is  shown  graphically. 
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TABLE  1 
SAMPLE  INPUT  DATA  SET  FOR  C/S  MODEL 


01/16/87  TEST  OF  (2X2),  1  UNSTABLE  ROOT  S 

NSTEP 

NTYPX  NTX1  NTX2  NTX3  loci    loc2   loc3  u{  u2  u3  r{  r2  r.3 

i\\i\22*\] 

d31cd32cd33 

cAl  CA2CA1> 
£21  ^22  ^23 
c31  c32c33 

NTYPY  NTY1  NTY2  NTY3  loci  loc2  loc3   v,  v0  v,  s,  s~,  s 


ab31ab32ab33 

a21  a22a23 
a31  a32a33 

DSEED  DSEED  DSEED 


1  v2  v3  »1  "2  >3 


TABLE  2 
SAMPLE  INPUT  DATA  SET  FOR  M/W  MODEL 


01/16/87  TEST  OF  (2X2),  1  UNSTABLE  ROOT  S 

NSTEP 

NTYPX  NTX1  NTX2  NTX3   Uj  u2  u3  r{  r2  r.3 

d11  d12d13 

d21  d22d23 
d31  d32d33 

£ll  £l2£l3 

£21  £22  £23 

c31  c32c33 

NTYPY  NTY1  NTY2  NTY3  w{  v2  v3  Sj  s2  s3 

fell  fel2&13 
521  P22  g23 
°31  ^32  P33 
all  a12a13 
a21  a22a23 
a31  a32a33 
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TABLE  3 
INPUT  DATA  SET  FOR  C/S  MODEL  CASE  1 

01/05/87  TEST  OF  (1X1),  NO  ROOTS 
99 

1  600  0  0  1.0  1.0  1.0  0.040  0.000  0.000  28  0  0 
0.06000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

1  400  0  0  1.0  1.0  1.0  0.010  0.000  0.000  16  0  0 
0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

345215789.D0      45635761. DO      89342761. DO 


TABLE  4 
INPUT  DATA  SET  FOR  M/W  MODEL  CASE  1 

01/05/87  TEST  OF  (1X1),  NO  ROOTS 
99 

1  600  0  0  1.0  1.0  1.0  0.040  0.000  0.000  28  0  0 
0.06000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

1  400  0  0  1.0  1.0  1.0  0.010  0.000  0.000  16  0  0 
0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
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Figure  E.l     Total  Force  Level  Trajectory  For  Case  One. 
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X  FORCE  LEVEL  TRAJECTORY 
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Figure  E.2    X  Force  Level  Trajectory  Over  Time  For  Case  One. 
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Figure  E.3    Y  Force  Level  Trajectory  Over  Time  For  Case  One. 
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TABLE  5 
INPUT  DATA  SET  FOR  C/S  MODEL  CASE  2 


01,05/87   TEST  OF  (1X1),  1  STABLE  ROOT  AT  200.200S 

99 
1  600  0  0  1.0  1.0  1.0  0.040  0.000  0.000  18  0  0 
0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

1  500  0  0  1.0  1.0  1.0  0.060  0.000  0.000  22  0  0 
0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761. DO   89342761. DO 


TABLE  6 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  2 

01/05/87  TEST  OF  (1X1),  1  STABLE  ROOT  AT  200,200S 
99 
1  600.     0.     0.   0.040  0.000  0.000     18      0     0 

0.00877  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00018  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

1  500.    0.     0.   0.060  0.000  0.000    22      0     0 

0.00874  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00017  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 
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TOTAL  FORCE  LEVEL  TRAJECTORY 
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Figure  E.4    Total  Force  Level  Trajectory  For  Case  Two. 
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X  FORCE  LEVEL  TRAJECTORY 
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Figure  E.5    X  Force  Level  Trajectory  Over  Time  For  Case  Two. 
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Y  FORCE  LEVEL  TRAJECTORY 
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Figure  E.6    Y  Force  Level  Trajectory  Over  Time  For  Case  Two. 
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TABLE  7 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  3 

01/06/87  TEST  OF  (1X1), 
99 
1    500      0     0    1.0    1.0 

1  UNSTABLE  ROOT  AT  200.200S 

1.0  0.010  0.000  0.000    22     0     0 

0.06000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

1    600      0      0    1.0    1.0 

1.0  0.010  0.000  0.000     18      0      0 

0.04000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

345215789. DO      45635761. DO      89342761. DO 

TABLE  8 
INPUT  DATA  SET  FOR  M/W  MODEL  CASE  3 

01/06/87  TEST  OF  (1X1),  1  UNSTABLE  ROOT  AT  200,200S 
99 

1  500.  0.  0.  0.010  0.000  0.000  22  0  0 
0.05129  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00018  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

1  600.  0.  0.  0.010  0.000  0.000  18  0  0 
0.03460  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00017  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
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TOTAL  FORCE  LEVEL  TRAJECTORY 
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TABLE  9 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  4 

01/05/87   TEST  OF  (1X1), 
99 
1    600     0      0    1.0    1.0 

ROOTS  AT 

"  200,200(U)  &  366,100(S)S 

1.0 

0.040 

0.000  0.000 

28 

0     0 

0.06000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

1    600      0      0    1.0    1.0 

1.0 

0.010 

0.000  0.000 

12 

0     0 

0.01000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

345215789. DO      45635761. DC 

•      89342761. DO 

TABLE  10 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  4 

01/05/87  TEST  OF  (1X1),  ROOTS  AT  200,200(U)  &  366,100(S)S 
1  600.     0.     0.   0.040  0.000  0.000    28      0     0 

0.04972  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00017  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

1  600.     0.     0.   0.010  0.000  0.000     12      0     0 

0.00796  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00017  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.10    Total  Force  Level  Trajectory  For  Case  Four. 
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TABLE  11 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  5 

01/31/87  TEST  OF  (1X2),  C 
99 
1    500  000     0    1.0    1.0 

ASE5 

1.0  0.030  0.000  0.000     38      0     0 

0.16000  0.06000  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00050  0.00040  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

2   250   300      0    1.0    1.0 

1.0  0.050  0.100  0.000    44     31      0 

0.08000  0.00000  0.00000 

0.06000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00000  0.00000 

0.00030  0.00000  0.00000 

0.00000  0.00000  0.00000 

345215789. DO      45635761. DO      89342761. DO 

TABLE  12 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  5 

02/02/87  TEST  OF  (1X2), 
99 
1  500.  000.     0.   0.030 

CASE  5 

0.000  0.000 

38 

0 

0 

0.12243  0.04733  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00038  0.00032  0.00000 

0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 

2  250.  300.     0.   0.050 

0.100  0.000 

44 

31 

0 

0.06156  0.00000  0.00000 

0.04530  0.00000  0.00000 

0.00000  0.00000  0.00000 

0.00015  0.00000  0.00000 

0.00023  0.00000  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.16    XI  Force  Level  Trajectory  Over  Time  For  Case  Five. 
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TABLE  13 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  6 

01/31/87  TEST  OF  (2X2),  C 
99 
2   200   300     0    1.0    1.0 

ASE  6,  NO  ROOTS 

1.0  0.010  0.020  0.000     11     16     0 

0.01000  0.02000  0.00000 

0.01000  0.01000  0.00000 

0.00000  0.00000  0.00000 

0.00008  0.00010  0.00000 

0.00010  0.00008  0.00000 

0.00000  0.00000  0.00000 

2   300  400      0    1.0    1.0 

1.0  0.005   0.010  0.000     15     15      0 

0.02000  0.02000  0.00000 

0.02000  0.01000  0.00000 

0.00000  0.00000  0.00000 

0.00010  0.00005  0.00000 

0.00005  0.00010  0.00000 

0.00000  0.00000  0.00000 

345215789.D0     45635761. DO      89342761. DO 

TABLE  14 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  6 

01/31/87  TEST  OF  (2X2), 
99 
2  200.  300.    0.   0.010 

CASE  6,  NO  ROOTS 

0.020  0.000 

11 

16 

0 

0.00935  0.01828  0.00000 

0.00902  0.00925  0.00000 

0.00000  0.00000  0.00000 

0.00007  0.00009  0.00000 

0.00009  0.00007  0.00000 

0.00000  0.00000  0.00000 

2  300.  400.    0.   0.005 

0.010  0.000 

15 

15 

0 

0.01671  0.01683  0.00000 

0.01638  0.00823  0.00000 

0.00000  0.00000  0.00000 

0.00009  0.00004  0.00000 

0.00004  0.00008  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.20    X  Force  Level  Trajectory  Over  Time  For  Case  Six. 
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Figure  E.21     Y  Force  Level  Trajectory  Over  Time  For  Case  Six. 
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Figure  E.22    XI  Force  Level  Trajectory  Over  Time  For  Case  Six. 
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Figure  E.23    X2  Force  Level  Trajectory  Over  Time  For  Case  Six. 
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TABLE  15 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  7 

01/31/87  TEST  OF  (2X2),  C 
99 
2   200   300     0    1.0    1.0 

ASE  7,  2  ROOTS 

1.0  0.010  0.020  0.000     11     16     0 

0.01000  0.02000  0.00000 

0.01000  0.01000  0.00000 

0.00000  0.00000  0.00000 

0.00008  0.00010  0.00000 

0.00010  0.00008  0.00000 

0.00000  0.00000  0.00000 

2   300  400      0    1.0    1.0 

1.0  0.005  0.010  0.000      9     12      0 

0.02000  0.02000  0.00000 

0.02000  0.01000  0.00000 

0.00000  0.00000  0.00000 

0.00010  0.00005  0.00000 

0.00005  0.00010  0.00000 

0.00000  0.00000  0.00000 

345215789.D0     45635761. DO      89342761. DO 

TABLE  16 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  7 

01/31/87  TEST  OF  (2X2), 

99 
2  200.  300.     0.   0.010 

CASE  7,  2  ROOTS 

0.020  0.000     11 

16 

0 

0.00944  0.01891  0.00000 

0.00942  0.00925  0.00000 

0.00000  0.00000  0.00000 

0.00007  0.00009  0.00000 

0.00009  0.00007  0.00000 

0.00000  0.00000  0.00000 

2  300.  400.     0.   0.005 

0.010  0.000      9 

12 

0 

0.01900  0.01905  0.00000 

0.01821  0.00957  0.00000 

0.00000  0.00000  0.00000 

0.00009  0.00005  0.00000 

0.00005  0.00009  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.26    Total  Force  Level  Trajectory  For  Case  Seven. 
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Figure  E.27    X  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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Figure  E.28    Y  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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Figure  E.29    XI  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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Figure  E.30    X2  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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Figure  E.31     Yl  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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Figure  E.32    Y2  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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TABLE  17 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  8 

02/01/87  TEST  OF  (2X2),  C 
99 
2   200   200     0    1.0    1.0 

ASE  8;  1  STABLE  ROOT 

1.0  0.100  0.120  0.000    51     85      0 

0.04000  0.06000  0.00000 

0.05000  0.07000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00030  0.00000 

0.00010  0.00040  0.00000 

0.00000  0.00000  0.00000 

2   300   300     0    1.0    1.0 

1.0  0.100  0.150  0.000    71     76     0 

0.06000  0.04000  0.00000 

0.05000  0.03000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00020  0.00000 

0.00030  0.00020  0.00000 

0.00000  0.00000  0.00000 

345215789.DO     45635761. DO      89342761. DO 

TABLE  18 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  8 

02/01/87  TEST  OF  (2X2), 
99 
2  200.  200.    0.   0.100 

CASE  8;  1  STABLE  ROOT 

0.120  0.000 

51 

85 

0 

0.03003  0.04312  0.00000 

0.03652  0.04845  0.00000 

0.00000  0.00000  0.00000 

0.00015  0.00021  0.00000 

0.00007  0.00028  0.00000 

0.00000  0.00000  0.00000 

2  300.  300.     0.   0.100 

0.150  0.000 

71 

76 

0 

0.04653  0.03115  0.00000 

0.03147  0.01856  0.00000 

0.00000  0.00000  0.00000 

0.00015  0.00016  0.00000 

0.00019  0.00013  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.33    Total  Force  Level  Trajectory  For  Case  Eight. 
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Figure  E.34    X  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  E.35    Y  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Figure  E.36    XI  Force  Level  Trajectory'  Over  Time  For  Case  Eight. 
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TABLE  19 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  9 

02/01/87  TEST  OF  (2X2),  C 
99 
2  400  400     0    1.0    1.0 

ASE  9,  1  UNSTABLE 

1.0  0.050  0.030  0.000    55    60     0 

0.08000  0.11000  0.00000 

0.12000  0.07000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00030  0.00000 

0.00020  0.00020  0.00000 

0.00000  0.00000  0.00000 

2   400  400      0    1.0    1.0 

1.0  0.070  0.020  0.000    65     54     0 

0.10000  0.05000  0.00000 

0.08000  0.15000  0.00000 

0.00000  0.00000  0.00000 

0.00020  0.00010  0.00000 

0.00010  0.00030  0.00000 

0.00000  0.00000  0.00000 

345215789. DO     45635761. DO      89342761. DO 

TABLE  20 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  9 

02/01/87   TEST  OF  (2X2), 
99 
2  400.  400.    0.   0.050 

CASE  9,  1  UNSTABLE  ROOT 

0.030  0.000 

55 

60 

0 

0.04814  0.06804  0.00000 

0.07278  0.04330  0.00000 

0.00000  0.00000  0.00000 

0.00012  0.00019  0.00000 

0.00012  0.00012  0.00000 

0.00000  0.00000  0.00000 

2  400.  400.    0.   0.070 

0.020  0.000 

65 

54 

0 

0.06706  0.03417  0.00000 

0.05404  0.10125  0.00000 

0.00000  0.00000  0.00000 

0.00013  0.00007  0.00000 

0.00007  0.00020  0.00000 

0.00000  0.00000  0.00000 
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Figure  E.40    Total  Force  Level  Trajectory  For  Case  Nine. 
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Figure  E.41     X  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Figure  E.42    Y  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Figure  E.43    XI  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Figure  E.44    X2  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Yl  FORCE  LEVEL  TRAJECTORY 
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Figure  E.45    Yl  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Figure  E.46    Y2  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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TABLE  21 

INPUT  DATA  SET  FOR  C,  S  MODEL  CASE  10 

02,02/87  TEST  OF  (3X3) 
99 
3   200    150    175    1.0    1.0 

1.0  0.005   0.010  0.020     4      4     11 

0.06000  0.02000  0.02000 

0.00000  0.00000  0.00000 

0.03000  0.01000  0.02000 

0.00000  0.00000  0.00000 

0.00000  0.00030  0.00020 

0.00000  0.00000  0.00010 

3    150    125   225    1.0    1.0 

1.0  0.003   0.020  0.030     4      5     14 

0.05000  0.01000  0.05000 

0.00000  0.00000  0.00000 

0.02000  0.01000  0.03000 

0.00000  0.00000  0.00000 

0.00000  0.00020  0.00030 

0.00000  0.00000  0.00010 

345215789. DO      45635761. DO      89342761. DO 

TABLE  22 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  10 

01/26/87   TEST  OF  (3X3) 
99 
3  200.  150.  175.   0.005 

0.010  0.020     4     4     11 

0.05836  0.01956  0.01916 

0.00000  0.00000  0.00000 

0.02701  0.00863  0.01745 

0.00000  0.00000  0.00000 

0.00000  0.00028  0.00019 

0.00000  0.00000  0.00009 

3  150.  125.  225.   0.003 

0.020  0.030     4      5     14 

0.04884  0.00993  0.04893 

0.00000  0.00000  0.00000 

0.01884  0.00926  0.02661 

0.00000  0.00000  0.00000 

0.00000  0.00019  0.00028 

0.00000  0.00000  0.00009 
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Figure  E.47    Total  Force  Level  Trajectory  For  Case  Ten. 
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Figure  E.49    Y  Force  Level  Trajectory  Over  Time  For  Case  Ten. 
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Figure  E.50    XI  Force  Level  Trajectory  Over  Time  For  Case  Ten. 


137 


X2  FORCE  LEVEL  TRAJECTORY 


o 
o 


o 


o 
o 


LEGEND 
C/S  MODEL 
M/W  MODEL 
C/S  PLUS  TS'ft" 
t/£  MINUS  1  STI 


T   I   I   II   1 1 1 1   I   II   I '  '  I 1   i   i   i 

S  10  15  20  25  30  35  40  45  50  55  60  03  70  70  80  65  00  93100 


TIME 


CASE  10  (3X3) 


Figure  E.51     X2  Force  Level  Trajectory  Over  Time  For  Case  Ten. 
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Figure  E.52    X3  Force  Level  Trajectory  Over  Time  For  Case  Ten. 
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Figure  E.53    Yl  Force  Level  Trajectory  Over  Time  For  Case  Ten. 
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Figure  E.54    Y2  Force  Level  Trajectory  Over  Time  For  Case  Ten. 
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